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Abstract 

Within the hmits of the large- A^c approximation (with the number of colors), we establish the 
high-energy behaviour of the diffractive and inclusive cross-sections for deep inelastic scattering 
at fixed impact parameter. We demonstrate that for sufficiently high energies and up to very large 
values of Q^, well above the proton average saturation momentum {Ql), the cross-sections are 
dominated by dense fluctuations in the target wavefunction, that is, by the relatively rare gluon 
configurations which are at saturation on the resolution scale of the virtual photon. This has 
important physical consequences, like the emergence of a new, diffusive, scaling, which replaces 
the 'geometric scaling' property characteristic of the mean field approximation. To establish this, 
we shall rely on a dipole version of the Good- Walker formula for diffraction (that we shall derive 
here in the context of DIS), together with the high-energy estimates for the dipole scattering 
amplitudes which follow from the recently established evolution equations with Pomeron loops 
and include the relevant fluctuations. We also find that, as a consequence of fluctuations, the 
diffractive cross-section at high energy is dominated by the elastic scattering of the quark- 
antiquark component of the virtual photon, up to relatively large virtualities ^ {QD- 
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1 Introduction 



In this paper, we shall present the first analysis of the physical consequences of the fluc- 
tuations in the high-energy evolution in QCD for the phenomenology of lepton-hadron 
deep inelastic scattering (DIS). By "fluctuations" we mean in general the correlations in 
the gluon distribution of an energetic hadron (the target in DIS) which are generated 
through the radiative processes contributing to its evolution with increasing energy. More 
precisely, we shall be interested here in the gluon-number fluctuations, which have been 
recently shown [1-4] to play an essential role in the evolution of the elastic S'-matrix for 
dipole-hadron scattering towards the unitarity limit. The "dipole" is a quark-antiquark 
pair in a colorless state, so hke the qq excitation through which the virtual photon couples 
to the hadronic target in DIS. As demonstrated in Refs. [2-4], both the energy depen- 
dence of the characteristic scale for the onset of unitarization — the target saturation 
momentum Qg — and the functional form of the scattering amplitudes at high energy, 
are strongly modified, and to a large extent even determined, by such fluctuations. Mo- 
tivated by these observations, and by the importance played by dipole configurations for 
the DIS processes at high energy, we shall provide a detailed analysis of the influence 
of gluon number fluctuations on the inclusive and diffractive DIS cross-sections. As we 
shall discover, the effects of fluctuations are in fact overwhelming: For sufficiently high 
energy, and within a wide kinematical range around the (average) saturation momentum, 
they wash out the corresponding predictions of the mean field approximation, and replace 
them with a new kind of universal behaviour, which is relatively simple. 

A priori, one may expect the gluon-number fiuctuations not to be important for the 
problem of the high-energy evolution, since the latter is characterized by a high-density 
environment (especially in the saturation regime at transverse momenta below Qs), where 
the gluon occupation numbers are large. And as a matter of facts, such fluctuations have 
been left out by the modern approaches to non-linear evolution in QCD, namely, the 
Balitsky equations for the eikonal scattering amplitudes [5] and the (functional) JIMWLK 
equation [6-8] for the evolution of the 'color glass condensate' [8,9]. These formahsms do 
include non-trivial correlations — as obvious from the fact that they generate infinite 
hierarchies of equations — , but these are purely color correlations (like the exchange of 
a single gluon between two dipole amplitudes) , and as such they die away when N^. oo 
(with Nf. the number of colors). In this limit, the Balitsky-JIMWLK hierarchy reduces to 
a single non-linear equation^ — the Bahtsky-Kovchegov (BK) equation [5, 10] — , which 
is often referred to as the mean field approximation (MFA). But in so far as the gluon- 
number fluctuations of interest here are concerned, the Balitsky-JIMWLK formalism as 
a whole is a kind of MFA — the generalization of the BK equation to finite Nc- 

Anticipated by early studies [11] (mostly numerical) within Mueller's "color dipole pic- 
ture" [12], the importance of the gluon-number fluctuations for the high-energy evolution 
in QCD has been fully appreciated only recently, over the last couple of years [1-4] . Rare 
fiuctuations associated with gluon splitting have been shown to control the formation of 

^ More precisely, they reduce to a simplified hierarchy which is equivalent to the BK equation 
for uncorrelated initial conditions, and which does not generate new correlations in the course 
of the evolution; see the discussion in Ref. [4]. 
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higher-point correlations in the dilute tail of the gluon distribution at relatively large 
transverse momenta (above Qs). Once initiated through fluctuations, such correlations 
are rapidly amplified by the BFKL evolution [13], and then play a crucial role in the 
growth of the gluon distribution and its eventual saturation: they are at the origin of the 
saturation and unitarization effects included in the Balitsky-JIMWLK equations. Thus, 
while they are correctly describing the role of the higher-point correlations in providing 
saturation, the Balitsky-JIMWLK equations fail to also include the physical source for 
such correlations, namely, the gluon-number fluctuations in the dilute regime. 

Soon after this failure has been flrst realized [4], new equations have been proposed [4, 14, 
15] which correct the Balitsky-JIMWLK hierarchy by including the proper source terms 
(the relevant fluctuations) within the limits of the large-iVg approximation. The latter 
is necessary since the dynamics leading to fluctuations in the dilute tail is described in 
the framework of the dipole picture [12, 16]. There is currently an intense, ongoing, effort 
aiming at generalizing these equations to flnite [17-24], but in spite of some interesting 
discoveries — like the recognition [18, 25] of a powerful self-duality of the high-energy 
evolution, and the construction of an effective Hamiltonian which is explicitly self-dual 
[20, 23] — , the general evolution equations valid for arbitrary Nc are not yet known. This 
is the main reason for our restriction to the large-A^c approximation in what follows. 

The new equations in Refs. [4, 14] form an inflnite hierarchy for the A'"-dipole amphtudes 
^j'(Ar)^ -^th N > 1 which describe the scattering between a system of N color dipoles 
and an energetic hadronic target. For reasons to be briefly mentioned in Sect. 2.1, we 
shall refer to these equations as the 'Pomeron loop equations'. They are currently under 
active investigation [25-32], and although their general solutions are not known, a lot 
of information is already available about their asymptotic behaviour at high energy, via 
the correspondence [3, 4, 33] between high-energy QCD and statistical physics (see the 
discussion below and also Sect. 3 for more details). This is the information that we shall 
use in our analysis of deep inelastic scattering throughout this paper. 

Previous studies of DIS at high energy, or small "Bjorken-x" [34-54] (see also the re- 
view papers [55-57]), have shown that diff'ractive DIS represents a privileged framework 
to study both gluon saturation — since the diffractive cross-section a^m appears to be 
dominated by relatively large qq configurations, with size r ~ i/Qs {Qs is the target satu- 
ration momentum, and increases with the energy as Ql{x) ~ — and fluctuations — 
since, as known since long [58, 59], the inelastic part of (Tdifr is a measure of the dispersion 
of the virtual photon wavefunction over its various Fock space components (like qq, qqg, 
etc.). Note, however, that the "fiuctuations" which appear in the standard discussion of 
diffraction [58, 59] refer to the wavefunction of the projectile — in DIS, the virtual photon 
and its various hadronic excitations — , and not to that of the target (the proton, or more 
generally some hadron h). In such previous studies, the target has always been treated in 
the spirit of the mean field approximation, like an optical potential or some fixed gluon dis- 
tribution, off which scatter the partonic components of the projectile. However, our main 
emphasis here will be precisely on the gluon-number fluctuations in the target wavefunc- 
tion, which are the relevant fluctuations for a study of the high-energy limit, and which 
influence all the DIS observables, and not only those associated with diffraction. Indeed, 
as we shall see in Sect. 2, the calculation of both inclusive and diffractive cross-sections 
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Fig. 1. "Phases" of the hadronic wavefunction in the kinematical plane Y-ln^Q"^ /Qq). When 

Y < Yds the ghion distribution scales with momentum in a 'geometric' way, i.e. it is a function of 
Q'^/Q^. When Y > Ids it scales in a 'diffusive' way, i.e. it becomes a function of ln{Q^ / Q^) / \/Y . 

in DIS at high energy and large naturally involves the amplitudes for dipole-target 
scattering, as determined by the Pomeron loop equations alluded to above. 

Let us clarify here what we shall mean by the 'high-energy regime of DIS' in what follows. 
This is the regime where, when increasing the energy (or decreasing x), the resolution scale 
of the virtual photon is simultaneously increased (although not as fast as the total 
energy, of course !), in such a way that remains comparable, or even larger, than the 
target saturation momentum Q^x). Indeed, the strict high-energy limit, in which the 
energy is increased at a fixed value for (and for a fixed impact parameter b), is less 
interesting, since conceptually clear: for sufficiently high energy, the cross-section reaches 
its unitarity (or 'black disk') limit, in which the projectile is completely absorbed by the 
target. (This behaviour will be manifest in the theoretical formalism that we shall use, 
which preserves unitarity.) On the other hand, the behaviour at relatively large Q^, such 
that > g2 

(x), is considerably more interesting, since in that region the scattering is 
weak, yet — as we shall discover — it is strongly influenced by the physics of saturation, 
within a large window in whose width is increasing with 1/x. Moreover, with increasing 
energy, the effects of fluctuations in the evolution of the target gluon distribution become 
more and more pronounced, and the interplay between fluctuations and saturation leads 
to a rich structure in DIS at high energy and ^ Qli.^)^ that we now explain. 

To that aim, let us refer to Fig. 1, which represents the kinematical plane Y- p (with 

Y = ln(l/x) and p = ln((5^/Qo)) where is a fixed reference scale) for DIS scattering, 
and also the physical phase-space for the QCD evolution of the hadronic target. Within 
the latter interpretation, Y is the target rapidity and can be assimilated to either 
the transverse momentum of a gluon inside the target {k± ~ Q), or the transverse size 
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of a dipolc within the projectile (r ~ l/<5)- The straight hne denoted as "Qf' in Fig. 1 
represents the average saturation hne in the target: {Ps)y = ^^i{Q^)/Qo) — -^^^ with A 
the 'saturation exponent'. This is an 'at'erag'e' since, in the presence of fluctuations, the 
saturation momentum becomes a random quantity: the stochastic evolution of the target 
gives rise to a statistical ensemble of gluon configurations, and the value of Ql can vary 
from one configuration to another. The dispersion a"^ = {pD — (Ps)^ which characterizes 
this variation will play a fundamental role in what follows. This stochastic process turns 
out to be a random walk in ps = In^Q'^/Ql), hence the dispersion rises linearly with Y 
[3] : cr^(y) ~ DirY, with D{r the 'front diffusion coefficient' (see Sect. 3). The second line 
in Fig. 1, denoted as "Q^", will be shortly explained. 

The most interesting region for us here is the scaling window in between these two lines. 
At lower values of Q^, to the left of the saturation hne, the gluons form a color glass con- 
densate with large occupation numbers ~ l/a^, the dipole amplitude reaches the unitarity 
limit (T) = 1, and the mean field description is well justified. At much larger Q^, on the 
right of the scaling region, one finds the standard perturbative regime, where the gluon 
system is very dilute, the dipole amplitude shows 'color transparency', (T(r)) oc r^, and 
neither saturation, nor fluctuations, do play any role. But in the scaling region at inter- 
mediate Q^, the gluon occupation numbers are relatively low and the scattering is weak, 
(T) <^ 1, yet the dynamics is generally influenced by both saturation and fluctuations, in 
proportions which depend upon the actual value of Y. 

The crucial scale for the present discussion is the rapidity Yds which separates between 
an intermediate energy regime at F < Ids ; where cr^ -C 1 and the mean-field picture 
roughly applies (with some important limitations, though, as we shall discuss in Sect. 3), 
and a high-energy regime at y > Ids ; where ^ 1 and the target expectation values 
are dominated by rare fluctuations up to very large values of p. (Parametric estimates for 
Ids and the other scales which appear in Fig. 1 will be given in Sect. 3.) 

Specifically, the lower part of the scaling region at F < Ids is the window for geometric 
scaling [33, 60-63]. In that region, the dipole amplitude and the gluon occupation number 
'scale' as functions of the dimensionless variable Q^/{Ql); e.g.. 



This behaviour is particularly interesting as it may explain a similar scaling observed in 
the HERA data for DIS at x < 0.01 and up to rather large (Q^ < 400 GeV^) [60]. 

But in the high-energy regime at F > Ids, the geometric scaling is progressively washed 
out by fluctuations [2] and eventually replaced (when Y ^ Ids) by a new type of scaling 
[3, 4], for which we propose the name of diffusive scaling: the dipole amplitude scales as a 
function of the dimensionless variable (p — {ps))/a. This scaling extends up to very large 
values of p, namely, it holds so long as p — (ps) <^ o"^. The borderline of this domain at 
PdiY) = (ps) + cr^ ~ (A -|- DiT-)Y represents the upper part of the line "Q^" in Fig. 1. The 
functional form of the dipole amplitude is known within this entire domain [4] (see Sect. 
3). In particular, deeply within the scaling window in Fig. 1, this is roughly a Gaussian: 



{T{r))y ^ (r'iQDr = exp[-7o(p - (p.»], with 70 



~ 0.63. 



(1.1) 





(1.2) 
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The physics behind this asymptotic form is remarkably simple and universal (i.e., insen- 
sitive to both the initial condition at low energy and the details of the evolution with 
increasing energy). It relies on just two basic facts: (i) saturation exists — meaning that, 
in the event-by-event description, T{p) = 1 for p smaller than some Ps, and (ii) the 
saturation scale Ps is a random variable distributed according to a Gaussian law, with 
expectation value (p^) and dispersion which both increase linearly with Y. Then, the 
scaling behaviour and the Gaussian shape manifest in Eq. (1.2) are a direct reflection 
of the corresponding properties of the probability distribution for pg. These properties 
transmit so directly from the probability law to the average amplitude since, within the 
range indicated in Eq. (1.2), the expectation value {T{p)) is small, {T{p)) <g 1, yet it is 
dominated by gluon configurations which are at saturation at the scale p of interest (i.e., 
for which Ps ^ p ^ (Ps), and thus T = 1) : Indeed, although relatively rare, such con- 
figurations yield contributions of order one, whereas the respective contributions of the 
^?/pzca/ configurations, for which ps ~ (ps), are exponentially suppressed (cf. Eq. (1.1)). 

We thus arrive at the important conclusion that, even though the scattering is weak on 
the average, meaning that the target looks typically dilute when probed on the scale r ~ 
l/Q, the average amplitude is nevertheless controlled by "black spots" within the target 
wavefunction, i.e., by rare fiuctuations with unusually large gluon density. At sufficiently 
high energy and up to extremely large values of (much larger than (Ql)), the physical 
quantities are dominated by saturated gluon configurations. This conclusion will be further 
substantiated by our analysis of the DIS cross-sections at high energy. 

Specifically, by using estimates like that in Eq. (1.2) for the dipole amplitudes (T^^^) 
(the general formulae of this type will be presented in Sect. 3), together with appropriate, 
'dipole', factorization fomualae for the inclusive and the diffractive DIS cross-sections 
(that we shall discuss in Sect. 2), we shall analytically compute in Sect. 4 the dominant 
behaviour of the respective cross-sections in the high-energy regime — by which we shall 
more precisely mean the diffusive scaling regime at ^ 1 and p — (ps) ^ o"^ (see Fig. 
1). The main conclusion of this analysis can be concisely formulated as follows: All the 
physical properties (hke the diffusive scahng or the Gaussian decrease in Eq. (1.2)) that 
have been previously noticed at the level of the dipole amplitudes translate literally to 
the DIS cross-sections. This is so since, as we shall demonstrate in Sect. 4, the various 
convolutions which relate the DIS cross-sections to the dipole amplitudes — namely, 
the convolutions with the BFKL kernel which enter the high-energy evolution of the 
projectile, and those with the wavefunction of the virtual photon — are dominated by 
small dipole sizes r ~ l/Q (within the kinematic regime of interest here). 

This last point is perhaps surprising, especially in relation with the diffractive sector, 
where the cross-section at large Q'^ {Q'^ » Q^) was known to be dominated by relatively 
large dipoles, with r ~ l/Qs [34,35,40,68]. This standard picture is correct indeed, but 
at high energy it applies only in the event-by-event description, and not also on the 
average. Indeed, even for 3> (Q^), the statistical ensemble representing the target 
does still contain gluon configurations ('black spots') for which ^ ; although rare, 
such configurations will dominate the (inclusive and diffractive) cross-sections, because 
the photon wavefunction strongly favors the small dipole sizes r ~ l/Q- 

Another surprise of our analysis is that, in the whole 'high-energy regime' defined as above 
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— i.e., up to rather large values of Q^, such that ^ {Ql) — , the diffractive cross- 
section is dominated by its elastic component, that is, by the elastic scattering between 
the onium as a whole and the hadronic target. (The 'onium' represents the wavefunction of 
the projectile at the time of scattering, as produced via the BFKL evolution of the original 
qq pair; see Sect. 2 for details.) Whereas expected in the 'black disk' regime at ^ (Qf) 
[56,58, 59]. this property may look surprising in the high-Q^ regime at ^ (Ql) (e.g., 
it contradicts the corresponding expectations of the MFA [40,47,50]), but it is in fact 
natural in view of our present wisdom: even for such high values of Q^, the cross-section 
is dominated by 'black spots' which are already in the 'black disk' regime. 

The dominance of the elastic scattering has interesting consequences for the dependence 
of the diffractive cross-section upon the rapidity gap Ygap : in the high-energy regime, 
the differential cross-section per unit rapidity gap is strongly peaked near l^ap = ^■ 
Alternatively — and this is the language that we shall prefer in this paper — the cross- 
section integrated over all the values of the rapidity gap Fgap from a minimal value Y^^^ 
up to Y is independent of the lower limit Y^^, so long as the latter is not too close to Y. 
This second language will be more convenient for us here since the approximations that we 
shall employ — a "leading logarithmic approximation" with respect to ln(l//3) = Y — Yg^p 

— will allow us to compute the integrated cross-section, but not also to control the details 
of the differential cross-section near P — 1. 

Let us mention here some limitations of our approach (besides the high-energy and large- 
Nc approximations already discussed). We shall not be able to describe the impact pa- 
rameter dependence of the scattering amplitudes; that is, we shall follow their evolu- 
tion with increasing y at a fixed impact parameter b. There are two reasons for that: 
First, the asymptotic solutions to the Pomeron loop equations are known only after a 
'coarse-graining' in impact parameter space, which eliminates the dependence upon b [4]. 
Second, more generally, we do not expect perturbation theory to properly describe the 
6-dependence of the scattering amplitudes, especially near the edge of the hadron, where 
the gluon density is low and the dynamics is sensitive to confinement. Because of that, in 
perturbation theory one can study the 'blackening' of the hadron disk, but not also the 
radial expansion of the 'black disk', i.e., the problem of the Froissart bound [64,65]. 

Furthermore, in our study of diffraction, we shall not be able to describe processes with 
arbitrary large diffractive mass Mx, or arbitrary small values of /? (see kinematics in Fig. 
2). Indeed, the formalism for diffraction that we shall develop in Sect. 2 can accommodate 
saturation effects in the wavefunction of the target (the hadron), but not also in that of 
the projectile (the virtual photon). This implies an upper limit on the rapidity ln(l//9) of 
the projectile, which is however parametrically large: ln(l//5) <^ (l/a^) ln(l/Q;^). Hence, 
in the case of diffraction, the 'high-energy hmit' that we shall consider is the limit in 
which the total rapidity Y and the rapidity gap Y^ap can be arbitrarily large, but such 
that their difference F — Fgap — ln(l//?) remains finite and constrained as indicated above. 

This paper is organized as follows: In the first half of it, which covers Sects. 2 and 3, we 
shall develop the formalism which is required for a calculation of the DIS cross-sections 
in the high-energy regime and for large N^, and in the second half (i.e.. Sect. 4), we shall 
explicitly compute these cross-sections, and compare to the respective mean-field results. 
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Specifically, in Sect. 2 we shall develop a factorization scheme for the diffractive dipole- 
hadron scattering at high energy (the 'dipole' being, of course, the qq pair produced 
via the dissociation of the virtual photon in DIS). Our final formula turns out to be an 
interesting combination between the 'color dipole picture' by Mueller [12] and an early 
formula for hadron-hadron diffraction by Good and Walker [58]: Within the frame in 
which the target rapidity coincides with the minimal rapidity gap ^g™™, the projectile 
is described as an onium — a collection of dipolcs produced via the BFKL evolution 
of the original qq pair over a rapidity interval Y — Y^^^ — , and the diffractive cross- 
section is expressed as the onium expectation value of the elastic cross-sections for all the 
dipole configurations within the onium. Whereas the emergence of such a formula is not 
unexpected (as the dipoles fulfill the main requirement of the analysis in Ref. [58], namely, 
they are eigenstates of the S'-matrix operator), our explicit derivation has the merit to 
clarify the importance of properly choosing the Lorentz frame in order for this formula 
to be valid. The relation to previous approaches [36,46, 50], in particular to a non-linear 
evolution equation proposed by Kovchegov and Levin [44], will be also discussed. 

The aforementioned factorization formula involves the average dipole amplitudes (T*^^^) 
evaluated at the target rapidity l^ap, which in our formalism can be arbitrarily large. In 
Sect. 3 we shall describe the calculation of these amplitudes in the high-energy regime, 
with the purpose of justifying expressions like Eq. (1.2) and completing the physical 
picture developed before, in relation with Fig. 1. The presentation in Sect. 3 is based on 
recent developments in Refs. [3, 14], but it also contains some new elements and conceptual 
clarifications, like a parametric estimate for Yds, the concept of diffusive scaling, and the 
recognition of an interesting 'rigidity' property for the average dipole amplitude in the 
presence of fluctuations. 

Sect. 4 is our main section, in which we establish the dominant behaviour of the DIS cross- 
sections in the high-energy regime. Remarkably, it turns out that this behaviour can be 
computed from the scattering amplitude (T(r))y for a single dipole — the original qq 
fluctuation of the virtual photon. This is a priori clear for the inclusive cross-section, but 
it is true for the diffractive one (at high energy), because of the dominance of the elastic 
scattering, as aforementioned. Our calculation will be organized as follows: First, we shall 
consider the qq component alone and compute the respective contributions to dtot and 
(Tdiff- Then, we consider the additional contributions associated with the qqg component 
and demonstrate that, in the regime of diffusive scaling, the inelastic contribution is 
parametrically suppressed as compared to the elastic one. To better emphasize the effects 
of ffuctuations, we shall compare our results to the corresponding predictions of the MFA. 

Sect. 5 summarizes our results and presents our conclusions and some perspectives. 



2 A dipole picture for deep inelastic scattering at high energy 

As explained in the Introduction, we would like to provide a theoretical description for 
diffractive deep inelastic scattering at high energy and in the multi-color limit ^ oo. 
By 'diffractive events' we shall understand the DIS events 7*/i — > Xp which contain an 
elastically scattered hadron h in their final state together with a 'diffractive hadronic 
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Y = ln(l/a;), 



gap — 



ln(l/a;p) 




ln(l//5) 



Fig. 2. Kinematics for diSractive DIS at high energy, or small Bjorken~x: is the virtuahty 
of 7*; W is the center~of-mass energy of the 7*/i system (with W"^ ^ Q"^)', M'j^ is the invariant 
mass squared of the diffractive system. 

state' X which is separated by a rapidity gap Ygap from the scattered hadron. (We use 
here standard notations in the context of DIS, as summarized in Fig. 2. See, e.g., Refs. 
[55, 56] for review papers on the phenomenology of diffraction and more details about the 
kinematics.) In the high energy regime at W"^ ^ Q^, M|-, the existence of a gap Ygap — 
ln[iy^/((5^ + M]^)] follows automatically from the condition that the proton undergoes 
elastic scattering together with the relevant kinematics. We shall be mainly interested here 
in the high-energy limit, defined as the limit W"^ ^ oo at fixed values for and Mj-. 
This means that the total rapidity difference Y = ln(l/x) between the projectile (7*) and 
the target (h), as well as the rapidity gap Ygap = ln(l/a;p) can increase arbitrarily large, 
but such that their difference Y — Ygap = ln(l//3) remains finite, and not too large (see 
Eq. (2.13) below for the precise condition). Note that our definitions for either diffraction 
or its high energy limit are not the most general possible ones (e.g., one could consider 
diffractive processes in which the target too breaks up after the collision), but they do 
cover interesting physical situations and, besides, they are constrained by the limitations 
of the subsequent theoretical developments. 

As also mentioned in the Introduction, we shall directly compute the diffractive cross- 
section integrated over the rapidity gap from some minimal value Yg™™ up to Y. For 
convenience, from now on, we shall reserve the simpler notation Ygap for this minimal 
rapidity gap, as this quantity will appear very often, and we shall measure the actual 
rapidity by the corresponding /5-variable, cf. Fig. 2. We have, clearly, the following relation 



between the integrated quantity in the l.h.s., that we shall explicitly compute in what 
follows, and the differential cross-section per unit rapidity — which appears under the 
integral sign in the r.h.s. and is the quantity usually considered in phenomenological 








(F,ln(l//?),Q2) 



(2.1) 
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studies of diffraction. As also indicated by the notations in Eq. (2.1), we shall consider 
cross-sections at fixed impact parameter b. The first, "elastic", term in the r.h.s. of the 
above equation is simply the boundary term for the integration at /5 = 1, or l^ap = ^ '■ 




{Y,Y,Q'), 



(2.2) 



and the reason why we refer to it as "elastic" is because, within the 'leading-logarithmic 
approximation' in ln(l//3) that we shall use, it corresponds to the elastic scattering be- 
tween the dissociation products of the virtual photon and the target (see Sect. 2.1 below 
for details and Appendix B for results going beyond the leading-logarithmic approxima- 
tion alluded to above). Note however that the scattering is not elastic from the viewpoint 
of the virtual photon itself (see the discussion at the end of Sect. 2.1 and also Fig. 8). 

Let us briefly describe here the physical picture underlying our subsequent developments. 
In a convenient 'dipole' frame in which the hadron h carries most of the total energy, the 
scattering between the virtual photon 7* and the target proceeds as follows: Long before 
the scattering (say, at time to — * ~oo), the virtual photon dissociates into a quark- 
antiquark pair in a color singlet state (a 'color dipole'), which then evolves through soft 
gluon radiation until it meets the hadron (at time t = 0) and scatters off the color flelds 
therein. At high energy, the color dipole and the accompanying soft gluons are eigenstates 
of the ^'-matrix operator — the collision acts merely as a color rotation on these states 
— , and the original picture of diffraction by Good and Walker [58] can be taken over: 
The diffractive process '~f*h Xh consists in the elastic scattering between the various 
Fock space components of the projectile and the target. More precisely, we shall argue 
below that this simple picture holds only in a well-tuned frame, in which the rapidity Yq 
of the target coincides with the (minimal) rapidity gap Y^ap in the final state ^ . Such a 
special choice of the frame is necessary in order to avoid the explicit treatment of the 
final state interactions, i.e., the gluon emissions and absorptions which take place in the 
wavefunction of the projectile after the time of scattering, and whose detailed description 
goes beyond the dipole picture. 

As implicit in the above picture, we shall restrict ourselves to the leading logarithmic 
approximation at high energy, in which the evolution consists in the emission of 'small- 
X gluons', that is, gluons which carry only a small fraction x <^ 1 oi the longitudinal 
momentum of their parent parton. So long as the energy is not too high, the evolution 
remains linear and is described by the BFKL formalism [13], which leads to a rapid growth 
of the gluon occupation numbers. But for sufficiently high energy, saturation effects which 
tame this growth start to be important [9,68,69], and introduce non-linear effects in 
the evolution equations [4-8, 10, 14]. In the physical situations that we shall consider, the 
projectile will be always dilute and thus evolve linearly. On the other hand, we shall allow 
for arbitrarily high energies (and thus for saturation effects) on the side of the target. 

The high-energy evolution will be further simplified through approximations valid at 
large Nc- On the projectile side, the combination of the BFKL evolution with the large- 



^ In practice, the identification of Yq with Ygap must hold within the accuracy of the leading 
logarithmic approximation: Yq and Y^ap can differ by an amount AY such that OLgdY <^ 1. 
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Nc approximation leads to Mueller's 'color dipole picture' [12] : When Nc oo, a gluon 
can be effectively replaced with a pointlike quark-ant iquark pair in a color octet state, 
and a soft gluon emission from a color dipole can be described as the splitting of the 
original dipole into two new dipoles with a common leg. In this picture, the original qq 
pair produced by the dissociation of the virtual photon evolves through successive dipole 
splittings and becomes an onium — i.e., a collection of dipoles — at the time of scattering. 
Strictly speaking, the dipole picture provides only the norm of the onium wavefunction, 
that is, the probability distribution for the dipole configurations within the onium [12, 16]. 
But, as we shall see, this knowledge is indeed sufficient to compute diffraction within the 
present formalism. On the target side, the large- A^c approximation turns out to be essential 
for the successful construction of evolution equations for the amplitudes describing the 
scattering between the projectile dipoles and the target gluons [4, 14, 15]. 



2. 1 Dipole factorization for the diffractive cross-section 

To develop our formalism, let us ignore for a while the electromagnetic process 7* qq (at 
high energy, this process factorizes out [70, 71] and can be easily reintroduced later; see the 
beginning of Sect. 4.2), and focus on the onium-hadron (Oh) scattering. The quantity 
that we would like to compute is the probability Pdiff(a;, ?/; y, Ygap) for the diffractive 
process Oh — > Xh in which the hadron h undergoes elastic scattering while the onium 
O dissociates into some arbitrary hadronic state X which is separated from the outgoing 
hadron by a rapidity gap equal to, or larger than, Fgap (with Fgap < Y, of course). In 
this definition, x and y are the transverse coordinates of the quark and, respectively, the 
antiquark in the original qq pair, the one which evolves into the onium. 

The probability P(iis{x,y;Y,Ygi,p) is, of course, frame-independent, but in what follows 
we shall derive an explicit expression for it by working in the specific frame in which 
Ygap coincides with the rapidity Yq of the target (and where the projectile has therefore a 
rapidity^ Y — Yq = Y — ^gap)- To avoid redundant notations, we shall replace everywhere 
Yga,p — >■ Yq and write, e.g., Pd[g{x, y; Y, Yq). The final formula that we shall arrive at (see 
Sect. 2.2 for a derivation, and Figs. 3 and 4 for graphical representations) reads 

Pdifr(a;,y;F,lo) = E dT^ PN{{z^}■,Y -Yq) {l - S{1)S{2) ■ ■ ■ S{N)) ,(2.3) 

N=l ^ ° 

with notations to be explained now: 

i) PN{{zi}]Y — Yq) is the probability density to produce a given configuration of N 
dipoles after a rapidity evolution Y — Yq starting with an original dipole {x,y). The 
configuration is specified by — 1 transverse coordinates {zj} = {zi, Z2, ...Zjv_i}, which 
physically represents the coordinates of the — 1 emitted gluons, and in terms of which 
the coordinates of the dipoles are {zq,Zi), {zi, Z2),---,{zn-i, Zn), with Zq = x and 



^ Within the leading-log approximation w.r.t. ln(l//3) = Y — Yq, we can ignore the difference 
in rapidity between the quark and the antiquark components of the original qq pair and treat 
that pair as a 'particle' with unambiguous rapidity Y — Yq. In Appendix B, we shall go beyond 
this leading-log approximation for the case where /? is close to one. 
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Fig. 3. Typical diagram contributing to the diffractive process j*h Xh in the frame where the 
target Yq coincides with the rapidity gap. For the projectile, we illustrate the gluon dynamics 
before and at the time of scattering. The final hadronic state X can be formed with an arbitrary 
number of gluons produced via 'final state interactions' (see the discussion in Sect. 2.2). The 
gluons in the target recombine back before the final state, so that the hadron emerges intact 
from the collision. For simplicity, we exhibit only two-gluon exchanges. 



Zjy = y. Also, J2n I dFAT with dT^ = d?Zid?Z2 . . .d^zj^^i represents the sum over all 
the configurations. The dipole probabilities are obtained by solving appropriate evolution 
equations (see below) with the following initial conditions : 

Pi(y = o) = i, p^>i(r = o) = o. (2.4) 

The evolution of the probabilities is such that the correct normalization condition 

oo „ 

^ / dTNPN{{zi}-^Y) = 1 (2.5) 

Af=l •' 

is satisfied at any Y . 

ii) S{i) = S{zi^i,Zi) is the S'-matrix for the scattering between the ith dipole in the 
projectile and a given configuration of color fields in the target. (Recall that, in the 
high-energy scattering, a dipole is an eigenstate of the interaction [70, 71] and the in- 
ternal configuration of the target is frozen during the duration of the collision; see also 
Eqs. (2.16)-(2.17) below.) We shall also need the scattering amplitude corresponding to 
a single dipole, defined as T{i) = 1 — S{i). Then: 

N N 

Sn = 1[S{i) and An = 1-Sn = l-Ul^-Tit)], (2.6) 

i=l i=l 

are respectively the S'-matrix and the scattering amplitude for a given configuration 
of N dipoles. The brackets in {An)yo denote the target average, that is, the average 
over the ensemble of color fields in the target. The target wavefunction, and thus the 
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Fig. 4. The same as in Fig. 3 hut in the large~Nc limit. Each gluon in the wavefunction of the 
virtual photon has been replaced with a pointlike quark-antiquark pair in a color octet state. The 
gluons inside the target wavefunction are not shown explicitly anymore. The relative simplicity 
of the qq representation allows us to also exhibit some multiple gluon exchanges, corresponding 
to unitarity corrections. 



corresponding expectation values, depend upon the rapidity interval Yq available for its 
internal evolution. Note that, in a given event (i.e., for a given configuration of the color 
fields in the target), the N dipoles scatter independently from each other — the total 
^-matrix iSat is simply the product of N factors corresponding to the individual dipoles 
— , but correlations are generally introduced by the average over the target (because the 
color fields there have non-trivial correlations). Thus, the expectation value {Sj\f)Yo is not 
factorizing anymore. Such target correlations will play an important role in the subsequent 
discussion in this paper. 

iii) The quantity |(1 — S{1)S(2) ■ ■ ■ S{N))yo\'^ = \{-^n)yo\'^ is recognized as the proba- 
bility for the elastic scattering between a given set of dipoles and the target. Thus, as 
anticipated, the diffractive probability (2.3) represents the projectile average of the elastic 
probabilities for all the possible (dipole) configurations in the projectile. Thus defined, 
the diffraction includes, but it does not reduces to, the truly elastic collision, in which 
the onium as a whole scatters elastically off the target. The difference appears because, 
whereas the individual dipoles (or gluons) are eigenstates of the interaction, this is not 
true also for their superposition (the onium), since the various states in this superposition 
interact differently with the target. 

The elastic probability Pel for onium-target scattering is rather computed as 

P,i{x,y;Y) = \1-S{x,y-Y)\' = \A{x,y-Y)\\ (2.7) 

where S{x,y]Y) is the diagonal ^-matrix element, S = {^in\S\'^in) = ml"^ out) , which 
measures the overlap between the final state emerging from the collision and the initial 
state prior to it. Therefore, Psurvivai = \<S\'^ is the probability that the original state (for 
the ensemble target+projectile) survives intact after the collision, whereas 

P,^,i{x,y;Y) = l-\S{x,y;Y)\\ (2.8) 
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is the probability for some inelastic process to occur. In what follows, the amplitude 
v4. = 1 — iS for the elastic scattering will be succinctly referred to as the forward amplitude. 
Within the present framework, this quantity can be computed as 

oo „ 

A{x, y;Y) = Y: / dr^ PN{{z^}; Y - Fq) (l - ^(1)^(2) ■ ■ ■ S{N)) , (2.9) 

which is perhaps a more familiar formula (various versions of it can be found in the 
literature [10, 30, 72]), and will be also derived below. 

As indicated by its notation, the forward amplitude A is independent of the rapidity 
divider Yq. This is a priori true on physical grounds — since, in the computation of A, 
Yq plays no dynamical role, but merely specifies the Lorentz frame — and is also verified 
by our explicit formula (2.9), within its accuracy limits (see the discussion below). In 
particular, one can choose to compute A in the frame where Yq ~ Y. In that frame, the 
projectile is just an elementary dipole {x,y), without additional gluons, and therefore 
A{x,y;Y) = 1 - {S{x,y))Y = {T{x,y))Y, in agreement with Eqs. (2.9) and (2.4). Simi- 
larly, the elastic probability (2.7) can be computed as the elastic scattering of the qq pair 
alone: according to Eqs. (2.3), (2.7), and (2.9) we can write 

Pei{x,y;Y) = \{T{x,y))Y\' = P^^{x,y-Y,Y) , (2.10) 

where the second equality reflects the physically obvious fact that an elastic scattering is 
the same as a diffractive event having Yq = Y. This explains the identification performed 
in Eq. (2.2). 

Via the optical theorem, the forward amplitude (2.9) also determines the total (or 'inclu- 
sive') cross-section^ : 

Ptot{x, y;Y) ^2 ReA{x, y; Y) = Pel + Pinei ■ (2.11) 

At this point, a word of caution is necessary, concerning a slight abuse in our termi- 
nology: The various "probabilities" introduced so far are strictly speaking differential 
cross-sections for onium-hadron scattering at fixed impact parameter; e.g., 

^(r,b,Y) = Ptot(x,y;Y), ^{r,b,Y,Yo) = PM^,y;Y,Yo), (2.12) 

(r = X ~ y and h = {x + y)/2 are the transverse size and the impact parameter of the 
original qq pair). As such, these quantities are certainly proportional to the corresponding 
scattering probabilities, but they are not necessarily bound to be smaller than one. Rather, 
they are constrained by the unitarity of the ^-matrix, which in the high energy regime 
(where the scattering amplitudes are predominantly real) requires < (An) < 1 for 
any A^. This condition, together with the above formulae, implies the standard inequality 
-Pdiff < |-Ptot [59], together with a series of upper bounds like ^ < 1, Pdifr < 1, and 



^ Note that our definitions for the scattering amplitudes T and A differ by a factor of i from the 
usual definitions in the textbooks. With our conventions, these amphtudes are predominantly 
real at high energy. 
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-Ptot < 2. Within the formaUsm that wc shall use below to compute the dipole amplitudes, 
all these constraints are correctly respected, and the various upper bounds are saturated 
in the high-energy limit. 

Before we turn to a derivation of Eqs. (2.3) and (2.9) in the next subsection, let us specify 
their range of validity in rapidity, and explain how to compute the dipole probabilities 
P/v(y — Yq) and the target-averaged matrix elements (^Ar)^;, which enter these formulae. 

The rapidity Y — Yq oi the projectile should be small enough for the saturation effects to 
remain negligible. This in turn requires [12, 16] 

1 N'^ 

Y-Yo^—\n^ (2.13) 
as al 

where dig = asNc/fi should be treated as a fixed quantity in the large-A^c limit. (Re- 
call that the typical rapidity interval necessary for the emission of one small-a; gluon is 
dY ~ l/as-) For rapidities satisfying this constraint, the probabilities Pn{Y — Yq) can be 
obtained by solving a 'Master equation' [16], actually, a set of coupled, linear, equations 
for the evolution with Y, whose structure makes it clear that the dipole evolution in the 
dilute regime is a Markovian stochastic process (see also Refs. [4,73]). Alternatively, the 
generating functional for P/V) to be introduced in Eq. (2.35) below, obey the non-linear 
evolution equation (2.36), originally derived by Mueller [12]. For latter convenience, let 
us display here the expressions for Pat generated after only one step in the evolution, that 
is, for y — Yq = dY and to linear order in the small quantity asdY : 



Pi{dY) ^1-dY^I dzM{x,y,z), 



P2{z\dY) = dY-^M{x,y,z), (2.14) 
Zn 

and Piv(d^) = for A?" > 3. In these equations, 

M{x, z) ^ , (2.15) 

[x - zy{y - zy 

is known as the 'dipole kernel' [12] : {as/'^'T^).Mix,y, z) is the differential probability for 
an elementary dipole {x,y) to split into two dipoles {x,z) and {z,y) per unit rapidity. 
Note that the integral over z in the formula for Pi has logarithmic singularities at z = x 
and z = y. Such singularities are expected at the level of the dipole probabilities, but they 
cancel out in the calculation of physical quantities, as it can be checked on the examples 
of Eqs. (2.3) and (2.9). 

Consider now the matrix elements for dipole-target scattering: At high energy, the 5*- 
matrix S{x, y) corresponding to a single dipole (a;, y) can be computed in the eikonal 
approximation as 

S{x,y) = j^iY{y\x)V{y)) , V\x) = Pexpjzy j cix-^+(x-, a;) i»| , (2.16) 

where the Wilson lines V\x) and V{y) represent the color rotations suffered by the quark 
and, respectively, the antiquark after their scattering off the color field in the target. 
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(The i"'s are the generators of the SU(A^c) algebra in the fundamental representation and 
the symbol P denotes path-ordering in .) Note that, in our conventions, the projectile 
propagates in the negative z (or positive x~) direction, so it couples to the /!+ component 
of the color field in the target. At high energy and in a suitable gauge, this is the only non- 
trivial component, and the average over the target wavefunction amounts to a functional 
average over (the 'color glass' average [8,9]) : 

{S{x,y))y, = / B[A+] WyM^] ^Jr(v\x)V{y)), (2.17) 

where the 'color glass weight function' iyyy[/l"'"] (a functional probability density) can be 
interpreted as the squared wavefunction of the target and depends upon the respective 
rapidity Yq. The expectation value {Sisr)Yo for the scattering of N dipoles is defined simi- 
larly. Note that, physically, the non-linear effects in included via the Wilson hnes in 
Eqs. (2.16)-(2.17) describe multiple scattering to all orders. 

In the high-gluon density regime where the gluon-number fluctuations become negligible, 
the evolution of the weight function Wyq with increasing Yq is described by the JIMWLK 
equation [6-8] — a functional, non-linear, equation of the Fokker-Plank type. Via equa- 
tions like (2.17), the JIMWLK equation generates an hierarchy of ordinary evolution 
equations for the dipole amplitudes {Sn)yo, originally derived by Balitsky [5]. The non- 
linear effects encoded in the JIMWLK equation describe gluon saturation in the target 
wavefunction and translate into unitarity corrections in the Balitsky equations. 

However, as pointed out in the Introduction, the gluon-number fluctuations missed by the 
Balitsky- JIMWLK equations are in fact essential for the physics at high energy, in that 
they represent the source for the higher-point correlations responsible for saturation. So 
far, the combined effects of fluctuations and saturation on the dipole scattering at high 
energy have been taken into account only in the large- A'c limit [2-4, 15], where a new set 
of evolution equations has been derived [4, 14] — the 'Pomeron loop equations' alluded to 
in the Introduction (see also Refs. [25,26,74]). In what follows, we shall assume that the 
dipole amplitudes {Sn)y() obey these new equations, which within the limits of the large- 
Nc approximation are valid up to arbitrarily high energy. In addition to the standard 
BFKL terms [13], these equations include non-/mear terms — which correspond to gluon 
saturation in the target and ensure the unitarization of the scattering amplitudes at high 
energy — and source terms — which correspond to dipole splitting in the dilute part of the 
target wavefunction and encode the relevant, gluon-number, fluctuations at large Nc- The 
perturbation theory for these equations can be organized in terms of "BFKL pomerons" 
(the Green's function for the BFKL equation [13]) which interact via "triple Pomeron 
vertices" [35, 74, 75] for Pomeron sphtting and merging. Accordingly, the solutions to 
these equations naturally encompass the Pomeron loops. In Fig. 5, we illustrate some 
Pomeron loops effects in the calculation of the forward scattering amphtude (2.9). 

To summarize, all the ingredients appearing in Eqs. (2.3) and (2.9) can be computed, at 
least in principle, by solving evolution equations which are explicitly known. These equa- 
tions are valid for large and up to rapidities Y and Yq which can be arbitrarily large, 
but such that their difference Y — Yq obeys the condition (2.13). Within this condition, 
the expression (2.9) for the forward amplitude is independent of the rapidity divider Yq, 
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Fig. 5. Pomeron loops in the forward amplitude for onium-hadron scattering. 

as it should. Indeed, it has been demonstrated in Refs. [30, 72] that the Fo^dependence in 
the r.h.s. of Eq. (2.9) cancels out exactly when using the Master equation for the dipole 
probabilities Pn{Y — Yo) together with the (large- A'^c version of the) Balitsky-JIMWLK 
equations for the dipole amplitudes {Sn)yo- Now, as argued before, the correct equa- 
tions at high energy are not the Balitsky-JIMWLK equations, but rather the Pomeron 
loop equations of Refs. [4, 14]. The latter include the effects of gluon number fluctuations, 
which are additional sources of lo^dependence. In a more general calculation of A{x, y; Y) 
which would be valid in any frame, such additional dependencies would be compensated 
by recombination effects in the wavefunction of the projectile. However, in any 'dipole 
frame' which satisfies the condition (2.13) such recombination effects are truly negligible, 
so Eq. (2.9) is indeed independent of Yq, up to higher-order corrections^. Incidentally, 
the above argument also shows that our formula for the diffractive probability, Eq. (2.3), 
is not independent of Yq, as expected on physical grounds. 

For later convenience, let us also introduce here the inelastic diffraction, a process which, 
as shown in Ref. [59] in the context of hadron-hadron collisions, is a direct probe of 
'parton' (here, dipole) fluctuations in the wavefunction of the projectile. The corresponding 
probability is obtained by simply subtracting out the (totally) elastic component (2.7) 
from the diffractive probability (2.3): 



P^^\x,y;Y,Yo) = PMx,y;Y,Yo) - P,,{x,y-Y) . 

This formula takes on a particularly suggestive form after noticing that Eqs. 
(2.9) can be rewritten as (with simplified notations) 



(2.18) 
(2.3) and 



^ This does not mean that the physical consequences of the gluon-number fluctuations are also 
neghgible (these will be discussed in the next sections). It is only the induced lo^dependence 
which is small indeed. 
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Pdiff(a;,2/;y,yo) = ( (^)targct )^^^., A(x,y;Y) ^ {{A)te.rget)^^^.. (2.19) 

In these formulae, the target-averaged amphtude ^target = (-A)targct is still an 'operator' 
from the point of view of the projectile, in the sense of depending upon a fixed configu- 
ration of dipoles. Also, at high energy, the amplitude is predominantly real, so one can 
ignore the modulus sign in the previous equations. We thus finally arrive at [59] 

P^i^: y, Fo) = ( A'a.get>p,„. ' { Aa.get)^^^. , (2.20) 

which clearly exhibits the fact that P^iff ' ^ measure of the dispersion of the dipole 
distribution within the wavefunction of the projectile. The dipole picture that we employ 
here provides an explicit realization for the projectile, and thus allows one to compute this 
dispersion and any other quantity pertinent to the distribution of dipoles. However, our 
main emphasis in what follows will be not on the effects of ffuctuations in the projectile 
(these are already well understood within the dipole picture; see, e.g., Refs. [4, 11, 12, 
16,72,73]), but rather on the physical consequences of the gluon number fluctuations 
in the target, as encoded in the evolution equations with Pomeron loops [4,14]. These 
fluctuations affect separately all the quantities introduced above {A, Pdm, -Ptot, etc.), 
because of their influence on the dipole scattering amplitudes {Sn)yo which enter the 
corresponding formulae. 

Let us conclude this subsection with a warning against the abusive interpretation of the 
notion of 'inelastic diffraction' in the context of deep inelastic scattering. As we shall see 
at the beginning of Sect. 4, the diffractive probabihty (2.3) determines the 'integrated' 
cross-section for DIS diffraction, i.e., the quantity in the l.h.s. of Eq. (2.1). In view of 
this, and of the identiflcation in Eq. (2.10), it becomes clear that the separation of adifr in 
between an 'elastic' plus an 'integral' piece, as shown in the r.h.s. of Eq. (2.1), corresponds 
to the decomposition of P^m in between its elastic and inelastic components, cf. Eq. (2.18): 
-Pdiff = -Pel + Pdm^- However, in DIS, and unhke in hadron-hadron coUisions, these two 
pieces of a^is cannot be separately measured. Indeed, even a process in which the onium 
scatters elastically (see Fig. 8) still appears as inelastic at the level of the DIS experiment, 
because the elastically scattered qq pair docs not recombine back into a virtual photon 
in the flnal state, but rather emerges as a hadronic state. Hence, the only experimentally 
relevant quantity is the 'total' diffractive cross-section, as determined by Pdiff. 



2.2 Justifying the dipole factorization 



To justify the above formulae for dipole factorization (in particular, those pertinent to 
diffraction) , we shall use the light-cone wavefunction formalism, as developed in the light- 
cone gauge A~ = (which is best suited for a study of the eikonal scattering of the 
projectile within our conventions). Our presentation will be rather sketchy, as similar 
techniques and manipulations can be widely found in the literature, and it will focus on 
the non-trivial aspects of the argument only. 

Let us denote by |\E'(— oo)), |\I'(0^)), |\I'(0"^)), and |^'(oo)) the wavefunction of the com- 
plete system target plus projectile (including their hadronic descendants in the case of 
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the final state) at times ^ t = — oo (long before the scattering), t = 0~ (just before the 
scattering), t = 0^ (immediately after the scattering), and, respectively, t = oo (the final 
hadronic state which is measured by the detector). We have: 

|*(0-)) = \0(Y - Yo)) \h(Yo)) , (2.21) 

where the two factors in the r.h.s. are the states of the projectile and, respectively, the 
target as produced after a rapidity evolution Y — Yq and, respectively, Yq. 

In the case of the projectile, the initial state for this evolution is the elementary color 
dipole {x,y) produced by the dissociation of the virtual photon: |C(0)) = \{x,y)) (the 
color indices are kept imphcit; see, e.g., Refs. [46,52] for details). The evolved state at 
t = 0~ is then a superposition of partonic states containing the original qq pair {x, y) 
together with an arbitrary number of soft gluons in a given spatial and color configuration. 
The detailed structure of such individual Fock states turns out not to be necessary for 
the present purposes. Rather, we shall simply denote by |A^) a generic state containing 
the qq pair {x, y) together with — 1 small-x gluons. We thus write, quite generically: 

\0{Y - Fo)) = E ^n{Y - Yo) \N) , (2.22) 

N 

where the sum over N should be really understood as a sum over the number of gluons, 
an integral over their transverse coordinates, and a sum over their polarizations and color 
configurations (we refer again to Refs. [46,52] for more precise notations). 

Concerning the target, its precise initial state at t = — oo is unimportant here. Rather, 
all that we need to assume is that at the time of scattering the target can be described as 
a superposition of states with a given color field (eigenstates of the gauge field operator 
A = A^) with coefficients $[A] which depend upon the rapidity Yq : 

\h{Yo)) = jD[A]^A]iYo) \A). (2.23) 

Consider now the onium-hadron collision which takes place at time t = 0. As anticipated, 
at high energy the partonic Fock-space states composing the onium are eigenstates of the 
collision operator: the transverse positions and the spins, or polarizations, of the 'partons' 
(quarks and gluons) are not changed by the scattering, while their color orientations un- 
dergo a field-dependent precession described by Wilson hues. We can write, schematically. 



1^(0+)) = / D[A] ^A]{Yo) E c^(F - Fo) (n S^]) \N) ^ \A) , (2.24) 

•' N i=l 

where Si[A] is the ^-matrix for the individual parton i, and is a Wilson line (with generally 
open color indices). 

^ Strictly speaking, the role of the 'time variables' in the present forniahsm is played by the 
hght-cone coordinates — x~ in the case of the projectile and, respectively, in the case of the 
target; for simplicity, we shall use the more intuitive notation t. 



20 



For all the processes that we are interested in, the hadronic target emerges intact from 
the collision, i.e., it undergoes elastic scattering, so we need only the projection of the 
outgoing state on the target state prior to the collision: 



l*diff(o+)) = |M>^o))(M^o)l*(o+)) 



N 



= / D[A] |$[A](yo)f E cn{Y-Yo) (n S^A]) \N) ® \h{Yo)) 

N i=l 

N 

^ E ^NiY - Fo) ( n S^[A]) \N) ® \hiYo)) , (2.25) 

N i=l ° 

where the brackets in the third line refer to the average over the target wavefunction. 
With the identification WVo[A] = |$[>l](yo)P, this target averaging is recognized as the 
color glass averaging in Eq. (2.17). 

Eq. (2.25) makes it manifest that, in the calculation of diffractive processes, the average 
over the target wavefunction is to be performed already at the level of the amplitude 
[39], rather than at the level of the squared amplitude (which is the quantity defining a 
probability; see Eq. (2.28) below). This peculiarity is, of course, related to our restriction to 
processes in which the target undergoes elastic scattering, and would not be true anymore 
for more general, diffractive, processes in which the target too is allowed to break up after 
the coUision (see also the discussion at the end of this subsection) . 

Note also that the target averaging in Eq. (2.25) automatically implies a color projection 
of the outgoing projectile onto color singlet states: Since the post-coUisional state of the 
target (which is the same as its initial state \h(YQ))) is a color singlet, so must be also the 
corresponding state of the projectile. In mathematical terms, the weight function l^y,JA] 
must be gauge-invariant, so the operation of averaging Oili Si[A] with this weight function 
must close all the Wilson hues into gauge-invariant color traces. Under this average and 
for large Nc, the product Dili Si[A] naturally reduces to the product of N dipolar factors, 
where is each of them is like the one appearing in Eq. (2.16). 

The final state at the time of detection is obtained by letting the 'diffractive' state in 
Eq. (2.25) evolve from i = up to i = oo under the action of the unitary evolution 
operator U (oo, 0) : 

|*diff(oo)) = C/(oo, 0) l^diff (0+)) . (2.26) 

It is this late time evolution which is responsible for the 'final state interactions' alluded to 
in the previous subsection: under the action of U (oo, 0), some of the soft gluons contained 
in the diffractive state at t — 0'^ may recombine with each other, so that the ensuing 
rapidity gap in a given event (that is, for a given final state) may be actually larger than 
Yq. (See Figs. 6 and 7 for an example.) In fact, the largest allowed gap is equal to Y and 
corresponds to elastic scattering^ , i.e., to the situation where all the gluons (or dipoles) 
in the final state recombine back before their detection into the original dipole {x,y). 
This is depicted in Fig. 8. 



This is consistent with the fact that, for Yq = Y, the diffractive probability (2.3) reduces to 

) inel 
diflf 



the corresponding elastic one, of. Eq. (2.10) : that is, P^^^{Y = Yq) = 0. 
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t=-co 00 -00 

Fig. 6. A particular diagram contributing to the difFractive process j*h Xh. The dashed verti- 
cal line at t = oo indicates the final (partonic) state. In the direct amplitude, a three-gluon state 
(together with the original qq pair) describes a component of the virtual photon wavefunction 
which interacts with the target hadron (not shown here) at time t = 0. After the scattering, final 
state interactions (gluon recombination) occur and at t = oo only two gluons are "measured". 
Similarly, in the complex conjugate amplitude, the wavefunction contains two soft gluons at the 
time of scattering; after the collision, these gluons recombine with each other, and finally an 
additional soft gluon is emitted before the time of measurement. The rapidity gap associated 
with this particular process is larger than Yq. 
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f = -oo 00 -00 

Fig. 7. The same as in Fig. 6 in the large-N^ limit; each gluon in the wavefunction of the virtual 
photon has been replaced with a pointlike quark-antiquark pair in a color octet state. 

On the other hand, within the present assumptions, the final gap can clearly not be smaller 
than Yq-. This is so because of our restriction to processes in which the target undergoes 
elastic scattering, so that none of the virtual quanta initially contained within the target 
wavefunction (and which are distributed in rapidity from to Yq) can be released in the 
final state. Since, on the other hand, there is no activity in the projectile wavefunction 
at rapidities smaller than Yq, it is clear that all the hadrons emerging in the final states 
and coming from the dissociation of the virtual photon must have a rapidity equal, or 
superior, to Yq. Therefore, the minimal rapidity gap is indeed Yq. (See Fig. 9.) 

In principle, it should be possible to compute the probability Pdig(^x,y;Y,Yg^p) for any 
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f = -oo 00 -00 

Fig. 8. Another diagram contributing to the difFractive process j*h Xh. Both the amph- 
tude and the c.c. amphtude correspond to contributions to the forward amplitude (2.9) for 
dipole-hadron scattering. Hence, this diagram is also an illustration of elastic dipole-hadron 
scattering, Eq. (2.7). (Notice that this is not an elastic scattering for 7* as well, since the virtual 
photon does not appear in the final state.) Here the rapidity gap takes its maximum value Y . 
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Fig. 9. A hnal example for the diffractive process j*h Xh. Here there are no final state 
interactions at all and the rapidity gap in this case is equal to its minimum value Yq. 



given value for the minimal rapidity gap Ygap within the range Y > Ygap > Yq. To that 
aim, one should explicitly follow the evolution of the outgoing state after the collision, 
as indicated in Eq. (2.26), and then project |\l/difT(c)o)) onto final states in which the 
hadrons released through the fragmentation of the projectile have rapidities Y > l^gap- 
This procedure would select those states within |\l/difr(0''')) in which all the projectile 
gluons with rapidities between Yq and Fgap have recombined back in the evolution from 
t = 0^ up to t = 00. Very recently, a more general formalism has been proposed [76], 
which in principle allows one to achieve this goal, and thus to compute the diffractive 
scattering in a generic frame (and for generic values of Nc). It would be interesting to see 
how our present factorization emerges from that formalism. 

However, our main point here is precisely that the calculation of P^if[{x, y; Y, Vgap) can be 
drastically simplified by conveniently choosing the Lorentz frame — namely, by choosing 
Yq = ^gap — , since with this particular choice one does not need to worry about the final 
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state interactions. Indeed, the condition that the gap be bigger than, or equal to, Yq is 
automatically satisfied within this kinematics, and it does not introduce any constraint 
on the final state of the projectile. Hence, the probability for such a diffractive event 
is measured by the norm (\l/difif(cxD)|\&diff(oo)) oi the final state, which is the same as the 
norm of the state |^difr(0''")) immediately after the collision, because of the unitarity of the 
operator f/(oo, 0). One sometimes says that 'the final state interactions cancel between the 
direct amplitude and the complex conjugate one', but in the present case this cancelation 
is an almost trivial consequence of our choice for the frame (see, e.g., Refs. [44, 77, 78] for 
less trivial cancelations of this type). 

More precisely, in computing P^is one needs to consider the deviation between the out- 
going diffractive state and the incoming state |\l/(0")) (indeed, the final state must be 
different from the initial one in order to have a real scattering). Thus, if one defines: 

|5^diff(0+)) = |^diff(0+)) - 1^(0-)), (2.27) 

then in this particular frame one can write 

Pam{x,y;Y,Yo) = (5^difr(0+)|^*d.fi(0+)) 

N 2 

-1:MY-Y,)\' {llS,[A]-l)^^ , (2.28) 



TV 



i=l 



where in writing the second line we have used the fact that the partonic states |A^) form 
an orthonormal basis. At large N^, the norm \c]sf{Y — lo)P of th^ state |A^) (summed over 
polarizations and color indices) can be identified with the probability Piv(^ — Yq) for a 
A^-dipole configuration in the dipole picture. With this identification, Eq. (2.28) is finally 
recognized as our originally proposed expression, Eq. (2.3). 

The expression (2.9) for the forward scattering amplitude can be similarly obtained : By 
definition, A = 1 — S with S the forward (S'-matrix element, computed as 

5(a;,2/;y) = (^(0-)|v^(0+)) 

= y: \cn{y - ro)r / D[A] \m{yo)\ n s^[A] , (2.29) 

N '' 1=1 

which immediately leads to Eq. (2.9), as anticipated. 

The previous discussion also explains our limitation to diffractive processes in which the 
hadronic target emerges intact from the collision. Of course, it would be very interesting 
(especially in view of applications to the phenomenology) to be able to describe the 
more general processes 7*/i — > XY in which the hadron dissociates into some low-mass 
hadronic system Y separated by a rapidity gap l^ap from the diffractive state X. But to 
that aim, one cannot avoid a detailed study of the final state interactions, which in the 
case of the target requires moreover a model for the hadronic structure. 

It is finally interesting to compare our formula (2.3) for the diffractive probability to 
related results in the literature. In the next subsection we shall show that, under an 
additional, mean field, approximation, our expression for Pdifr obeys a non-linear equation 
originally proposed by Kovchegov and Levin [44]. But before that we would like to point 
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out that Eq. (2.3) also encompasses other results from the previous literature. When 
the target is not too dense so that multiple scattering is negligible, Eq. (2.3) reduces to 
the expression employed by Bialas and Peschanski [36]. Namely, in the single scattering 
approximation, we have An ~ J2iLiT{i) (c.f., Eq. (2.6)), and then Eq. (2.3) becomes 

PMx,y;Y,Yo)^ I n{u,v;Y-Yo) \{T{u,v))yo\'' 

+ I n(')(ui,Vi;w2,'y2;l^-y'o)(T(wi,-yi))y„(r(u2,V2))^„,(2.30) 

Ui,Vi 

where n and n^^^ denote, respectively, the average dipole number density and the aver- 
age density of (distinct) dipole pairs in the projectile (see, e.g.. Sect. 5 Ref. [4] for the 
corresponding definitions) and the amplitude {T)yq obeys the BFKL equation, as appro- 
priate for a dilute target. Hence, Eq. (2.30) represents the BFKL approximation to Pdiff- 
Physically, the first term in the r.h.s. of Eq. (2.30) represents the probability that the 
same dipole scatters in both the direct and the complex conjugate amplitude, whereas 
the second term describes the scattering of different dipoles. 

On the other hand, for a generic, dense, target but a relatively small rapidity for the 
projectile, such that (5^(1^ — Iq) 1) one can restrict oneself to onium configurations 
which contain only two dipoles at the time of scattering. In that case, Eq. (2.3) reduces 
to results previously obtained by Kovner and Wiedemann [46] and by Munier and Shoshi 
[50]. The corresponding expressions will be presented in Sect. 4.2, where their high-energy 
limit will be also investigated. 

Let us emphasize here, however, some important differences in our treatment of the tar- 
get expectation values as compared to Refs. [46, 50]. (i) In Ref. [50], the target averages 
have been estimated in a mean field approximation (cf. Sect. 2.3) which neglects the cor- 
relations induced by gluon-number fluctuations; such an approximation is reasonable at 
intermediate energies, but it eventually fails at sufficiently high energies (see the discus- 
sion in Sect. 3). (ii) In an attempt to relax the restriction to elastic scattering on the 
target side, the authors of Ref. [46] have suggested to replace the target averaging at the 
level of the amplitude with an averaging at the level of the (diffractive) probability. With 
our present notations, their suggestion amounts to the following replacement 

(^)lof {APsin,A*)Yo (2.31) 

within the r.h.s. of Eq. (2.3). Here, Pging is the projector onto color singlet states for the 
projectile wavefunction, and it has been introduced to ensure that A involves only gauge- 
invariant operators; e.g., Pging-^ = An, cf. Eq. (2.6), for a A^-dipole state. (In our previous 
developments, such a projection was automatically ensured by the target averaging at the 
level of the amplitude, cf. Eq. (2.25).) However, with this new prescription, the target 
averaging is tantamount to summing over all the possible gauge-invariant final states for 
the target, without any restriction on their distribution in rapidity. But prior to scattering, 
the target wavefunction had developed virtual excitations at all the rapidities y ranging 
from y = up to y = Yq, and in the absence of any rapidity veto on the final state, 
there is no reason why these quanta should not materialize into hadrons occupying this 
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whole rapidity interval. In other terms, with the target expectations values evaluated as 
in Ref. [46], one is actually including processes which have no rapidity gap whatsoever. 
Such processes do not qualify as "diffractive" according to the usual terminology, and will 
be not considered in what follows. 



2.3 Evolution equations in the mean field approximation 



Although our main focus in this paper will be on the effects of fluctuations (i.e., of the 
deviations from the mean field behaviour) in so far as the target expectation values are 
concerned, it is nevertheless interesting at this level to slightly deviate from the main 
stream of the presentation and discuss the mean-field version of our precedent results. 
This is useful for, at least, two reasons: First, it will allow us to make contact with a 
non-linear equation for the diffractive probability previously proposed by Kovchegov and 
Levin [44], and thus clarify the conditions for the validity of the latter. Second, the mean 
field approximation (MFA) that we shall introduce here will later serve us as a term of 
comparison, to better emphasize the consequences of the fluctuations at high energy. 

Specifically, the MFA consists in the following factorization assumption 

{S {1)3(2) . . . S{N))y, ^ {S{1))yo {S{2))y, ■ ■ ■ (S(iV))y„ , (2.32) 

which neglects the correlations among the dipoles induced by their scattering off the 
target. Ultimately, this is an assumption about the absence of correlations in the target 
gluon distribution. (For instance, such a factorization holds indeed, at large N^., within 
the framework of the McLerran-Venugopalan model [9] for the gluon distribution of a 
large nucleus. More generally, it amounts to use a Gaussian approximation for the CGC 
weight function in Eq. (2.17) [79].) In what follows we shall demonstrate that, with this 
additional assumption, our formula (2.3) for the diffractive probability obeys indeed to the 
Kovchegov-Levin (KL) equation of Ref. [44]. During this procedure, and under similar 
assumptions, we shall also provide a rapid derivation of the Balitsky-Kovchegov (BK) 
equation for the forward amplitude [5, 10]. 

For more clarity, we shall denote with a bar quantities computed in the mean field ap- 
proximation. For instance: 

- °° r 2 

PMx,y;Y,Yo) = E dr^ PN{{zi};Y -Yo) l-s,S2---Sn , (2.33) 

N=l '' 

and similarly 

oo „ 

S{x, y;Y) = / dr^, PN{{zi}; Y - Fq) ^1^2 • • • Sn • (2.34) 

AT 1 



N=l 



In these formulae, Si is a simplified notation for {S{zi-i, Zi))Yo, which in the present context 
is of course the same as S{zi_i, z,; yo)- 

In what follows, we shall exploit the evolution of the probabilities Pn{Y — Yq) within the 
dipole picture in order to deduce a set of non-linear evolution equations for the quantities 
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defined in Eqs. (2.33)-(2.34). To that aim, it is more convenient to use the original version 
of the equations for P/v, due to Mueller [12]. Specifically, Mueller has derived an equation 
for the following generating functional 



oo „ 

Z^y[Y,u] = dr^ Pn{zi...zn-i]Y)uiU2 - ■ -UN, (2.35) 

N=l •' 

where Ui = u{Zi_i,Zi) is an arbitrary 'source' field, and {x,y) = {zq,Zn)-, as usual. The 
dipole probabilities Pn{X) can be deduced from Za,y[y, li] by functionally differentiating 
with respect to u and then letting — > 0. The evolution equation satisfied by Zxy reads 



dZxy 
dY 



M{x, y, z) i-Z^y + Z^^Z^y) , (2.36) 



to be solved with the following initial condition, which follows from Eq. (2.4) : 

Zxy[0,u] = Uxy. (2.37) 

As demonstrated in Ref. [73], this non-linear equation is equivalent with the infinite 
hierarchy of linear, master, equations for the probability densities P/v [16], although the 
respective pictures of the evolution are quite different: According to Eq. (2.36), the one- 
step evolution consists in the splitting of the original, high-rapidity, dipole into two child 
dipoles, which then separately evolve and produce their own distribution of dipoles. By 
contrast, the master equations of Ref. [16] focus on the splitting of the low-rapidity dipoles 
produced in the previous steps of the evolution. 

Eq. (2.36) looks formally similar to the BK equation for the -S'-matrix [10], and indeed 

the latter can be easily derived from the former, as we show now: Eqs. (2.34) and (2.35) 
imply S{x,y;Y) = Zxy[Y,s], which together with Eq. (2.36) immediately leads to (with 
the simplified notation Sxy = S{x, y; Y)) 



dSxy 

dY 



^jMix, y, z) {-Sxy + SxzS^y) , (2.38) 



which is the BK equation [5, 10], as anticipated. It is more customary to write this equation 
in terms of the scattering amplitude Txy = 1 — Sxy, in which case it reads: 



dTxy Otg 

dY ~ 2^ 



j M{x, y, z) (fx, + f,y - fxy - fxzf.y) . (2.39) 



The evolution of the diffractive probability (2.33) can be similarly addressed. Eqs. (2.33) 
and (2.35) imply (the dependence upon the rapidity gap Yq is kept implicit, since this 
variable is fixed in the following manipulations) 

P^^{Y) = l-2Sxy{Y) + Zxy%s'] 

= -1 + 2fxy{Y) + Zxy[Y, s% (2.40) 
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where we have assumed that s is real, as appropriate at high energy. This rewriting of 
P^y together with the previous equations (2.36) and (2.39) then imply 



z 

jM{x,y,z) [Ziyis'] - ZUs']Z,y[s']) 

z 

M{x, y, z) (2f^, + 2T,y - 2T^y - 2f,,f, 



2nJ 



zy 



-{P^^ + 1 - 2f^y) + {P^f + 1 - 2f^^){P^^ + 1 - 2f^y) ) (2.41) 



After some simple manipulations, the expression in the r.h.s. can be recast into the form 



dPi^ 



dY ^ i / ^' {^^^ + - + ^^^^ 

z 

-2f^,P^^ - 2P^^f,y + 2f ^,f,y) , (2.42) 



which is recognized, as anticipated, as the equation proposed by Kovchegov and Levin 
[44]. It is interesting to notice that, in Ref. [44], this equation has been obtained by 
working in the target rest frame, which required a more intricate analysis of the final 
state interactions. 

The non — linear terms in Eqs. (2.39) and (2.42) describe (incoherent) multiple scattering 
between the dipoles in the projectile and the color fields in the target and are responsible 
for unit arizat ion: In the high energy limit, the solutions T{Y) and Pdiff(^)^o) approach 
the 'black-disk' fixed points T — P^ss = 1, in agreement with the corresponding properties 
of the more general formulae (2.3) and (2.9). 

But although they do respect the unitarity bounds, Eqs. (2.39) and (2.42) cannot be 
used in a study of the high-energy limit, because of the mean field approximation (2.32) 
inherent in their derivation and which fails at high energy [1, 2]. To properly include the 
relevant fluctuations, one must replace the BK equation with the hierarchy of Pomeron 
loop equations [4, 14] for the iV-dipole amplitudes {S{1)S{2) ■ ■ ■ S{N))y- Once these am- 
plitudes are thus computed, they can be used to evaluate the diffractive probability ac- 
cording to Eq. (2.3), which replaces the solution to the KL equation for sufficiently high 
energy. Although the general solution to the Pomeron loop equations is not known, the 
correspondence between high-energy QCD and problems in statistical physics [3, 33] has 
allowed one to deduce valuable information about the behaviour of the dipole amplitudes 
in the high-energy limit [2-4] . This will be explained in the next section. 
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3 Dipole amplitudes at high energy: Fluctuations & Diffusive scaling 

The essential ingredient required by a calculation of the various scattering probabilities 
introduced in the previous section are the A^-dipole amplitudes {T^^^)y which describe the 
scattering between the target and a system of dipoles. In this section, we shall describe the 
calculation of these amplitudes in the high-energy regime where the MFA breaks down, 
because of the strong influence of gluon-number fluctuations [2-4] . 

Most of the results to be presented below have already appeared in Refs. [2-4], but 
since these are recent developments and, moreover, are of utmost importance for the 
present analysis, our respective discussion here will be quite complete. In doing so, we 
shall also clarify some points which have not been addressed in the previous studies, 
like the borderline between the intermediate-energy regime, where the BK equation is 
a reasonable approximation and geometric scaling applies, and the high-energy regime, 
where the evolution is dominated by fluctuations leading to diffusive scaling. 

The subsequent picture will be given in a frame in which all of the energy is carried 
by the hadronic target, so the projectile is a bare dipole, or a set of few such dipoles. 
To compute a target-averaged amplitude like (T'^^))y, one needs to (i) evaluate the 
scattering amplitude Ty(a;,y) for a single dipole and in a single event (meaning, for a 
given configuration of the gluon fields in the target, as generated through a single evolution 
from 1/ = up to 1/ = F) and (ii) average the result over the ensemble of events (that 
is, over all the possible target evolutions from y — to y — Y). It turns out that, for 
sufficiently large values of Y , both these operations become relatively simple and their 
results are universal (i.e., independent of the initial conditions at low energy), as they are 
fully determined by the quantum evolution with Y [2-4]. In fact, this universality is even 
stronger, in the sense that the details of the evolution matter for a couple of parameters 
which enter the final results (namely, the coefficients A and Dfr which appeared in the 
Introduction) , but not for the functional form of the average amplitudes [4] . 

3.1 The event-by-event dipole amplitude 

From now on, we shall neglect any non-trivial dynamics in the impact parameter space, 
that is, we shall assume that the evolution with increasing energy is quasi-local in b and we 
shall often omit the 6-argument from the amplitudes. This is a reasonable approximation 
so long as one is interested (as we are here) in the high-energy limit of the DIS cross- 
sections at fixed impact parameter. Then, the dipole amplitude in a single event can be 
written as 

TY{x,y) = Ty(r) = 7V(p) (3.1) 

where r = |£C — is the dipole size and p represents r in logarithmic units, like in Fig. 
1 : p = ln(l/r^QQ), with Qq a scale of reference introduced by the initial conditions at 
low energy. Note that large p corresponds to small dipole sizes, or to large transverse 
momenta (p ~ In /c^) after a Fourier transform. 
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For sufficiently fiigfi energy, tfie amplitude takes the form of a traveling wave [33] : 

Ty{p) ^T{p-p,{Y)). (3.2) 



This is a front which interpolates between a strong scattering regime at p < PsiX)-, where 
the unitarity bound T = 1 has been saturated, and a weak scattering regime at p > PsiX)-, 
where Ty <^ 1, and which propagates towards larger values of p when increasing Y . Here, 
Ps {Y) is the position of the front — conventionally defined as the dipole size p for which 
T[p) — 1/2 — and it increases linearly with Y: ps{y) — \oisY . The position of the front 
also serves as a definition for the target saturation momentum : 

Ps{Y) ^ \ri{Ql{Y)/Qt) ^ Ql{Y) ^ Q^e^-^. (3.3) 

Accordingly, A will be alternatively referred to as the saturation exponent^ or the front 
velocity. The BK equation (2.39) predicts A = Aq ~ 4.88 [61,62], but this result is signif- 
icantly lowered by gluon number fluctuations [2, 3] (see below) and also by the next-to- 
leading order corrections in perturbative QCD [63] . 

We shall also need later some analytic control over the shape of an individual front. As 
indicated in Eq. (3.2), this depends only upon the difference p — Ps(Y), meaning that 
the front propagates without distortion: in a co-moving frame, the shape of the front 
is independent of Y . This property has been first noticed within the framework of the 
mean field approximation, that is, for the solution to the BK equation [33,61,62,80-84], 
but for the individual fronts this remains true also in the presence of fiuctuations [3,66]. 
Since p — Ps — ln(l/r^(5^), we deduce that the dipole amplitude Ty(r) depends upon 
its two kinematical variables r and Y only trough the dimcnsionless product r'^Ql^Y) — 
the property usually referred to as geometric scaling [33,60-63]. This property has far 
reaching physical consequences, as it provides a natural explanation [61] for a new scaling 
behaviour observed in the small-x data at HERA [60] and plays also an essential role in 
understanding the particle production in deuteron-gold collisions at RHIC [85,86]. Note, 
however, that the quantity which matters for the calculation of physical observables is 
the target- averaged dipole amplitude {T{p))y and, as we shall later explain, the geomet- 
ric scaling property characteristic of a single event does not necessarily translate to the 
average quantities, because of fluctuations [2-4]. 

The fiuctuations are inherent in the evolution with increasing energy, because the instan- 
taneous configuration of the gluons in the target (as probed by the external dipole) can 
abruptly change from one rapidity step to another via the emission of new gluons. At 
high energy, the fluctuations are relatively important only in the dilute regime at large 
transverse momenta {k_[_ 3> Qs)) where the gluon occupation numbers are small. Since the 
high-A;_|_ gluons scatter preponderantly with small dipoles, of size r ~ 1/A;_l, we deduce 
that the fiuctuations should significantly affect only the tail of the dipole amplitude Ty (p) 
at large values of p. Simple physical considerations [3, 66] (see also below) show that the 
fluctuations have a relatively small effect on the overall shape of an individual front — 
they modify only the foremost part of the front, where T is very small: T ~ — , but 
they strongly influence the dynamics of the front, thus considerably reducing its velocity 
A with respect to the mean field prediction Aq. This behaviour, which is confirmed by 
numerical simulations within both QCD and statistical physics (see, e.g., [28,29,66,87]), 
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reflects the pulled mLtnie of the saturation front: the propagation is entirely driven by the 
dynamics in the tail of the front at p ^ PsO^), and thus is very sensitive to even small 
changes in this tail due to fluctuations. 

Returning to a discussion of the shape of the front, it is important to stress here (in 
anticipation of Sect. 3.2) that, in order to compute the average amplitudes at high energy, 
one needs only a very limited information about the shape of the individual fronts: it 
is in fact enough to know that 7V(p) = 1 for p < Ps{Y) [4]. But precisely in order to 
develop the corresponding argument in Sect. 3.2, and also in view of numerical tests in 
the forthcoming sections, we shall need to know a little bit more about the function 
Ty(p). Specifically, the following piecewise approximation, which can be inferred from the 
relevant literature [3,66,67], will be sufficient for our purposes: 

1 for 2; < 

Ae-'^°' for l<z<L (3.4) 
B e'"" for L <^ z. 

Here, z = p — Ps is the scaling variable, A and B are undetermined normalization factors, 
7o ~ 0.63, and L is the 'width of the front', namely, the distance p — Ps over which the 
amplitude falls off from its saturation value T = 1 to a value of order a^, where the 
fiuctuations become important. This condition immediately yields: 

L ~ - In ^ + 0(1). (3.5) 
70 aj 

Note that the physical regions in which the solution is approximately given have no overlap 
with each other, so, not surprisingly, the expressions shown in Eq. (3.4) do not provide 
a continuous interpolation for T{z). The first line in Eq. (3.4) is, clearly, the saturation 
region ; the second line describes the front of the traveling wave, with the characteristic 
'anomalous dimension' 1 — 70 ~ 0.37 which is a hallmark of the BFKL evolution in the 
presence of saturation [61, 62]; finally, the third line represents the perturbative QCD tail, 
where the amplitude exhibits color transparency (T(r) ~ r^), and which is generated via 
the large-p behaviour of the BFKL kernel. More precisely, within the BFKL evolution, 
the color transparency is amended by quantum corrections in the "double-logarithmic 
approximation" (DLA). These are not explicitly shown here, since the contribution of 
this large-p tail to physical quantities is anyway negligible in the high-energy regime of 
interest. It is however important to notice that such DLA effects violate geometric scaling, 
which is therefore restricted to the region at z ^ L. 

The single-event amplitude in Eq. (3.4) looks quite similar to the solution to the BK 
equation (see, e.g., Ref. [86]), but it differs from the latter in two essential points: 

a) The front is compact in Eq. (3.4), i.e., it has a finite width L which is independent of Y. 

By contrast, the front of the BK solution extends up to distances p — Ps ^ VY, and thus 
its width increases with Y. In the context of the BK equation, this front is often referred 
to as the geometric scaling window [61,62]. We see that, as a consequence of fiuctuations, 
the width of the geometric scaling window for a single event (and for sufficiently high 



T{z) = 



31 



energy) is considerably reduced with respect to the mean field approximation. 

b) The front velocity A which enters Eq. (3.4) via Ps is considerably smaller than the 
corresponding BK velocity Aq. This is seen numerically for generic values of a^, and is 
confirmed by an analytic calculation valid in the weak coupling limit a1 0, which shows 
that in this limit A converges towards Aq, but only very slowly [2,3,66] : 

C 

- >'o - 1 2/-, / 2^ < 1 . (3.6) 

In [l/aj) 

Here, C is a known number, which is determined by the linear, BFKL, dynamics and turns 
out to be quite huge : C = 7r^7ox"(7o) ~ 150. (Here, ^(7) is the Mellin transform of the 
BFKL kernel [2,3]. Note also that Aq = x(7o)/7o ~ 4.88.) The tendency of A to decrease 
when increasing is also confirmed by a study of the 'strong noise limit', which shows 
that in that hmit A vanishes as a power of l/a^ [88]. 

More precisely, the above picture applies only for high enough rapidities, larger than the 
formation time of the front, estimated as [3, 89] : 

as inform ~ 7r-;T7 — (3.7) 
2x" (70)70 

This is the rapidity evolution necessary for the front to reach its final form in Eq. (3.4), 
and for the velocity to reach its asymptotic value in Eq. (3.6), starting with some generic 
initial condition at F = 0. For Y < Fform, the width of the front increases with Y as 
P - Ps ^ via BFKL diffusion, while for Y > Iform, this width gets stuck at its 

maximal value, equal to L. As indicated by the 'geometric scaling' region in the diagram 
in Fig. 1, a similar behaviour is shown also by the average amplitude {T{p))y, but only 
up to some maximal rapidity Y ~ Ids, beyond which the effects of fluctuations become 
overwhelming (see the discussion in the next subsection). 

Let us conclude this presentation with a few physical considerations on the role of fluc- 
tuations [3]. The previous arguments exhibit the special role played by the parameter 
in the study of fluctuations (in particular, the mean fleld limit is obtained as — > 0), 
and this deserves a comment. In the dilute regime where the fluctuations are important, 
and for large Nc, the target itself can be described as a collection of dipoles, and then 
the amplitude T(r) describes the scattering between the external dipolc r and the target 
dipoles. The dipolc-dipole amplitude T(r, r') is of order and is peaked at sizes r' ~ r. 
Hence, to a good approximation, T(r) ~ alf{r), where /(r) is the dipole occupation num- 
ber in the target, that is, the number of dipoles with size r (per unit of In 1/r^) within an 
area around the impact parameter of the external dipole. In a given event, / is discrete 
(/ = 0,1,2,...), so the scattering amplitude takes on only discrete values, which are 
multiples of a^. In particular, the minimal non-zero value for T is a^, showing that the 
tail of the front must abruptly end when T becomes of O(a^). This explains the compact 
nature of the front. 

Furthermore, the particle number fluctuations follow a normal distribution: 6f ~ ^/f, 



showing that 6T ~ tt^vT ~ yO^jT. Thus, the mean fleld approximation becomes reliable, 
in the sense that ST <^ T, as soon as T ^ a^. In particular, this becomes a good 
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approximation everywhere (i.e., for any r) in the hmit — > 0. 



3.2 Front dispersion through fluctuations 



Let us consider now the average over the target wavefunction, which in the present context 
amounts to an average over the statistical ensemble of fronts : the fronts associated to all 
the possible evolutions of the target over a rapidity interval Y. It turns out that at high 
energy — sufficiently high for the fronts to reach their canonical form, cf. Eq. (3.7) — this 
averaging is quite simple: Since all the fronts in the ensemble have the same shape, they 
can differ from each other only by a translation. That is, the only random variable in the 
problem is the position ps of the front, which can be argued to be distributed according 
to the following, Gaussian, probability density [3, 66, 67] 



Py{Ps) = exp 
Vvrcr 



{ps - {ps)f 

(72 



(3.8) 



In this equation, (p^) ~ XagY is the average position of the front and increases with Y 
with the velocity shown in Eq. (3.6), since this is the common average velocity for all the 
fronts in the ensemble. Furthermore, = {pD — (ps)^ is the front dispersion, which 
rises linearly with Y, cr2(F) ~ DfragY, since the stochastic process is a random walk 
around the average front ^° . It has been first suggested by numerical simulations [66] , and 
very recently confirmed through analytic arguments [67] that, in the weak coupling limit 
al ^ 0, the front diffusion coefficient Df,. scales as 

In (1/q;2) 

This vanishes, as expected, in the mean field limit a'^ 0, but only very slowly. The 
coefficient V has been explicitly computed in Ref. [67]. (Once again, this is determined 
by the linear dynamics in the tail.) 

Then the average amplitude (T)y is determined by 



{T{p))y = I dp, Py{Ps) T{p - p,), (3.10) 

—oo 

with T(p — Ps) the single-event front in Eq. (3.4). Higher-point correlations can be com- 
puted similarly. Note that the ensuing average amplitudes will naturally depend upon the 
difference z = p — (ps), but they will also show additional dependencies upon Y, via the 
front dispersion a. That is, unlike the individual fronts, the average amplitudes will gen- 
erally not show geometric scahng [2, 3]. In fact, as we shall shortly discover, at sufficiently 
high energies the geometric scaling is replaced by a new, diffusive, scaling. 



" Note that our present normalizations for the coefficients A and Df^ differ by a factor a, from 
those used in the Introduction. 



33 



Given the Gaussian nature of the probabihty distribution (3.8), it is straightforward to 
compute the average in Eq. (3.10) for the piecewise-defined front in Eq. (3.4), and we 
shall do so indeed in the more detailed studies in the next subsections. Here, however, 
we would like to point out (following an original analysis in Ref. [4]) that, for sufficiently 
high energy, the dominant contributions to the average amplitudes {T^^^)y with any N 
are given by the saturation region in Eq. (3.4), and thus they are all expressed through a 
universal function. 

Let us start by better specifying what we mean by "sufficiently high energy" : As antici- 
pated in the Introduction (cf. Fig. 1), this is the regime in which the target expectation 
values are dominated by dense fluctuations, which in turn requires relatively high values 
of Y and not too large values of p (but p can still be much larger than the average sat- 
uration scale (ps); see below). More precisely, we shall demonstrate here that within the 
high-energy regime defined as 

High-energy regime <^=^ a ^ I/70 and p — (p^) <^ 7ocr^ , (3-11) 

the dominant contributions to the average dipole amplitudes can be simply computed 
with a ©-function saturation front: 



Tg{p) = e{ps-p), (3.12) 
which in turn implies [4] 

{T{p))y = ^Erfc(^) , (3.13) 
and, more generally, 

{T{p,mp,)...T{p^))y = ^Erfc(^) , (3.14) 

where z — p — (ps) and 2;> = p> — (p^), with p> = the largest among the p^'s. In these 
equations, Erfc(x) is the complementary error function. 



Erfc(a;) = ^ [ dte"*', (3.15) 
^/^^ J 

X 

which arises here as the integral of the probability density (3.8) over all values of Ps larger 
than p (respectively, p>): The r.h.s. of Eq. (3.13) is simply the probability to find a front 
which is at saturation on the resolution scale p. 

Eqs. (3.13)-(3.14) show that the average scattering amplitudes at high energy are domi- 
nated by those fronts within the statistical ensemble at Y which are at saturation at the 
highest resolution scale p> set by the incoming projectile. This is true, in particular, for 
relatively large values of p, well above (ps), where on the average the scattering is weak, 
{T{p))y <^ 1, yet this average value is dominated by the relatively rare configurations 
('black spots') which are at saturation for that (large) value of p. The subdominant con- 
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tributions not shown in Eqs. (3.13)-(3.14) are suppressed by, at least, one power of l/cx 
and/or z/a^ (respectively, z>/o"^), and are sensitive to the detailed shape of the single- 
event front. A similar sensitivity holds for the average amplitudes outside the high-energy 
regime defined in Eq. (3.11). 

Note that the universal behaviour exhibited in Eqs. (3.13)-(3.14) emerges also in the limit 
of strong fluctuations, as recently shown in Ref . [88] . This observation is consistent with 
the analysis in Ref. [4] and the present considerations, since the limit of strong fluctuations 
is formally obtained by letting o" ^ oo within our formulae, in which case the 'high-energy 
regime' of Eq. (3.11) extends everywhere. 

To demonstrate the above results, we need to study the influence of the shape of the front 
on the average amplitudes with increasing energy. To that aim, it is sufficient to consider 
the following, simplified, profile for an individual event : 



T,{p) = e(p, - p) + e(p - ps) e-^^'^-'^^^ . 



(3.16) 



For 7 = 7o a; 0.63, this simulates the front region in the second line of Eq. (3.4), while 
for 7 = 1 it rather corresponds to the 'color transparency' region in the third line there. 
By using Eqs. (3.16) and (3.8) it is straightforward to deduce 



m^)>y = ^Erfc(^)+iexp(^ 



7^2 



2-Erfc(^-^ 
a 2 



(3.17) 



with z = p— (ps). The ffist term is the same as in Eq. (3.13) and comes from the saturation 
piece of T^, while the second term comes from the exponential piece of Eq. (3.16). 

The following hmiting behaviours of the function Erfc(x) will be useful in what follows: 



Erfcfx) 



2 - 
1 

exp 



exp(— x^) 

y/7T\x\ 



for X <^ —1 
for X — 
for X ^ 1. 



(3.18) 



i) Consider ffist the (relatively) low energy situation at cr <^ I/7. By using Eqs. (3.17)- 
(3.18) one finds that the average amplitude retains its single event profile, except in the 
short interval |2;| ^ cr where it gets smoothed: 



(Tiz)) 



Y 



1 for 2; ^ — c 
i"'^^ for a <^ z, 



for cr <^ 1/7 . 



(3.19) 



In this regime, geometric scaling is manifest at the level of the average amplitudes. The 
situation is similar to the mean field approximation, except for the facts that the real front 
is compact and has a smaller velocity, cf. Eq. (3.6). (The compactness becomes manifest. 



Such estimates should more properly read, e.g., but here and in what follows we shall 

often use the fact that 70 ~ 0.63 is a number of 0(1) to simplify the parametric estimates. 
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of course, only if one starts with the more reahstic profile in Eq. (3.4).) We thus conclude 



In^fl/a^) 

Negligible dispersion <(=^ cr -C I/70 agY < ^' ^' . (3.20) 

The upper limit above, namely cksIds ~ is parametrically larger than the 

'critical' rapidity a^Yc ~ 1ii(1/cKs) ^'^'^ onset of the unitarity corrections, cf. Eq. (2.13), 
and also than the formation time agYforni ~ ^"^{^/(^l)^ ■'^Q- (3-7), for the individual 
fronts (as it should for the present discussion to make sense). 

ii) For a 3> I/7, the first term in Eq. (3.17) dominates everywhere except at extremely 
large distances ahead of the front, such that z ^ 70"^. For instance, in the interesting 
range cr -C 2; -C 7cr^, where the average amplitude is small, we obtain 

{T{z))y ~ - |l + ^1 e-^'/-' for a « ^ « (3.21) 
zy TT z y ^(T J 

where the first term inside the curly brackets comes from the 'saturation' piece in the 
r.h.s. of Eq. (3.17), while the second term corresponds to the 'exponential' piece there, 
and is suppressed with respect to the first term by a factor zj^o"^ ^ 1, as anticipated. 
One can similarly check that Eq. (3.14) yields indeed the dominant behaviour for the 
A^-dipole amplitude (T*^^))y within the range specified by Eq. (3.11). 

On the other hand, it is easy to check that, for much larger values of p, such that z ^ 
cr^, the dominant contribution to the average amplitude comes from the second term in 
Eq. (3.17), that is, from the exponential tail in Eq. (3.16). One then finds 

{T{z))y ~ exp(-7^) exp (^^) ^ > ^^ (3-22) 

where we recognize the same exponential decay with p as for a single front, Eq. (3.16). 
The piecewise expression of the average amplitude at high energy reads therefore 

'(l/2)Erfc(;2/(T) for -00 < z -C 7oO-^ 



Y 



for a > 1/70, (3.23) 

~ for 7ocr^ ^ 



where we have directly considered the more realistic profile in Eq. (3.4), so the exponential 
tail at very large z is determined by the 'color transparency' piece of Eq. (3.4). 

For all but very large values of 2;, the average front profile in Eq. (3.23) has no resemblance 
at all with the single-event profile in Eq. (3.4): There is no trace of the BFKL 'anomalous 
dimension' 70, nor of 'geometric scaling'. Rather, so long as 2; <^ 7oC^; the average ampli- 
tude obeys a new type of scaling — it scales as a function oi zja — which was previously 
observed in Refs. [3, 4] , and for which we propose the name of diffusive scaling . Note 
that the kinematical region for diffusive scaling coincides with the high-energy regime of 
Eq. (3.11), and is represented as the region at F > Ids and p at the left of the curve 
"Qd" in Fig. 1. In particular, within the weak scattering regime at cr <^ 2; <^ 7ocr^, the 
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Fig. 10. The average dipole-hadron forward amplitude of Eq. (3.17) (with 7 = Ij as a function 
of the "standard" scahng variable r^{Ql) and for various values of the front dispersion a. 
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Fig. 11. The same as in Fig. 10 as a function of the "diffusive" scaling variable 
z = — \n(r'^ (QD) / a . As a increases the amplitude approaches the limiting behavior in Eq. (3.13). 



average amplitude decreases as a Gaussian in z = p — (ps), and not as an exponential (as 
it happens for the single-event front within the geometric scaling window). 

This behaviour is illustrated in Figs. 10 and 11, which show the average amplitude in 
Eq. (3.17) for 7 = 1 and various values of cr, as a function of either the 'geometric scaling' 
variable r'^{Ql) (in Fig. 10), or the 'diffusive scaling' variable z = \r(1 / {Ql)) / a (in Fig. 
11). The "mean field" curve in Fig. 10 represents the hmit cr — > 0, that is, the single-event 
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profile in Eq. (3.16). This shows geometric scahng, and would be the actual behaviour in 
the absence of fluctuations. As manifest on Figs. 10 and 11, when increasing c, the average 
amplitude deviates more and more from geometric scaling and eventually approaches the 
hmiting profile in Eq. (3.13) (displayed as the "Erfc" curve in Fig. 11), which shows 
diffusive scaling. 

These considerations show that the mean field approximations (like the BK equation) 
completely fail to describe the dipole scattering at high energy. This is particularly striking 
when one considers multi-dipole correlations : in the high-energy regime (3.11) these 
correlations are given by Eq. (3.14), which shows that the mean-fleld factorization (2.32) 
is maximally violated. For instance, within the range cr -C 2; -C cr^, the average amplitudes 
are weak, yet strongly correlated : 

{Tip)) ^ (T'ip)) ^ (t^(p)) ^ ^ ^ e-^^/'^^ for a « ^ « (3.24) 

Accordingly, the A'^-dipolc amplitudes with enter the evolution equations with Pomeron 
loops become as important as the one-dipole amplitude (T(p))y already in the weak 
scattering regime, and not only in the region where the unitarity corrections are important. 
This is so because the individual fronts which contribute to these average amplitudes are 
themselves at saturation, and thus they are sensitive to unitarity corrections in the event- 
by-event description. 

This dominance of saturated gluon configuration within the diffusive scaling region has 
another interesting consequence, to which we shall refer as the rigidity of the average 
amplitude: This is the property that, within the intermediate weak scattering regime at 
a <^ z a^, the average amplitude increases unusually slow when increasing Y at fixed 
dipole size (i.e., at fixed p). To better appreciate this property, recall first the rapid, BFKL, 
increase of the average amplitude in the weak scattering regime at low and intermediate 
energies, say, in the geometric scaling window (cf. Eqs. (1-1)): 



d{T{r))y 



~ -foXas{T{r))y when r^Q^,) < 1. (3.25) 



dY 

On the other hand, Eq. (3.13) implies 

d{T{r)) 
dY 




Aa^ ^ (T(r))^ , (3.26) 



where we have used z = p — (ps), {ps) — Aa^F and ~ D^^agY, and the last estimate 
holds for (7 <^ z <C cr^ (compare to Eq. (3.24)). Clearly, as compared to the intermediate 
energy (or mean-field) regime in Eq. (3.25), the relative variation of the high-energy 
amplitude in Eq. (3.26) is smaller by a factor z/a^ <^ 1. This is so since the individual 
fronts which contribute to {T{r))Y in this regime are all at saturation, so by themselves 
they cannot increase when further increasing Y; the comparatively slow growth visible in 
Eq. (3.26) is rather due to the progression of the average front towards larger values of p. 
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4 Fluctuation effects on deep inelastic scattering at high energy 



In this section, we shall consider the high-energy limit of the inclusive and diffractive 
cross-sections introduced in Sect. 2, with the purpose of demonstrating the fundamental 
change of behaviour introduced in this regime by the gluon-number fluctuations. On this 
occasion, we shall present the first calculation of the DIS cross-sections including the 
effects of the Pomeron loops, and we shall compare our results with the corresponding 
predictions of the mean field approximation (the BK equation) . 



4.1 The general set-up 



To perform this analysis, we need to complete our previous factorization formulae with 
the electromagnetic vertex describing the dissociation of the virtual photon into a qq pair. 
Within the present, high-energy, approximations, this vertex factorizcs out and the DIS 
cross-sections for the (inclusive and diffractive) 7*/i scattering can be computed as [70, 71] 



da? 



tot 



{Y,Q')^ £dvjd'r Y \raiv,r;Q)\'Ptotib,r;Y), (4.1) 



d% ^ ^ TT 

a=L,l 



and, respectively (recall that Yg^p denotes the minima/ rapidity gap), 

{Y,Y,,^,Q')^ ['dv [d'r Yl mv,r;Q)\'Pais{b,r;Y,Y,,^). (4.2) 



d% 



a=L,T 

In these equations, are probability densities for the qq dissociation of a virtual 

photon with cither transverse (T) or longitudinal (L) polarization (these can be found 
in the literature and arc displayed, for convenience, in Appendix A), r = x — y and 
b — {x + y)/2 are the transverse size and the impact parameter of the qq pair, r = \r\, 
V is the fraction of the photon longitudinal momentum taken away by the quark, and 
Ptot and Pdiflf are the total and, respectively, diffractive probabilities for onium-hadron 
scattering, as introduced in Eqs. (2.3) and (2.11). 

The above formula are a priori frame-independent, but in order to be able to use our 
expression (2.3) for Pdiffj we need to evaluate the diffractive cross-section (4.2) in the 
frame where Yq = Vgap- On the other hand, the inclusive cross-section is most simply 
evaluated in the frame where Yq — Y and the projectile is an elementary qq pair: 



dotot 
d% 



f dv [ dV Y Q) r 2 Re {T(b, r))y . (4.3) 



a=L,T 



The factorization formula (4.2) holds within the same approximations as previously used 
in the calculation of Pdifr, namely, the leading logarithmic approximations with respect to 
both ln(l/a;p) = Ygap and ln(l//3) = Y — ^gap- Under these assumptions, the diffractive 
cross-section depends upon xp and (3 only via the high-energy evolutions of the target and 
the projectile, respectively. In the high-energy regime of interest, these approximations 
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are certainly correct in so far as the xp-dependence is concerned (since xp is truly small), 
but they fail to describe the actual /3-dependence in the vicinity of = 1 (i.e., in the 
case of a relatively small mass ^ for the diffractive system). In fact, within the 
present approximations, the diffractive cross-section for /3 ~ 1 reduces to its 'elastic' piece 
in Eqs. (2.1) and (2.2), which is independent of /3. However, for /3 ~ 1, the projectile is 
simply a qq pair, and for that pair the generalization of Eq. (4.2) which provides the correct 
/9-dependence is known [36,37] and will be exhibited in Appendix B (see Eq. (B.l)). It 
is therefore quite straightforward to improve the subsequent analysis in the vicinity of 
/? = 1, whenever this is needed. However, since our main interest here is rather in the 
high-energy limit — > at fixed /3, we find it more convenient to focus on generic values 
of /3 <^ 1, for which the relatively simple factorization formula (4.2) apphes as written. 

By using Pdiff = -Pel + P^^-> cf- Eq. (2.18), it is further possible to separate the 'elastic' 
and 'inelastic' components of the diffractive cross-section (4.2), as shown in Eq. (2.1). 
Although these components cannot be separately measured in DIS, as explained at the 
end of Sect. 2.1, this decomposition is still useful because, as we shall argue below, the 
elastic com'ponent dominates at high energy, up to relatively large values of (within 
the region for diffusive scaling) . Namely, we shall find that within the high-energy region 
defined by Eq. (3.11), the diffractive probability is predominantly elastic: 

Pdiff (r; F, Fgap) ~ Pei(r; Y) when (770 > 1 and In — ^ < a"^. (4.4) 

Since, moreover, the convolutions with the virtual photon wavefunction in Eqs. (4.1) and 
(4.2) are dominated by dipole sizes r ~ l/Q within the whole kinematical region of interest 
(this will be explicitly checked below) , the dominance of the elastic scattering extends to 
the diffractive cross-section, as anticipated : 

^(l^,nap,Q') ^ when a7o»l and In ^ « (4.5) 

This in turn implies that, within the present approximations at high energy, o"difr is quasi- 
independent of the minimal rapidity gap Ygap and can be simply computed as the elastic 
scattering of the elementary qq pair in the frame where Yq = Y (cf. Eq. (2.10)): 

drr^ /•! r 



f dv fd'r ^ mv,r;Q)\' \{T{b,r))y\\ (4.6) 



In view of this and of Eq. (4.3), we conclude that, for sufficiently high energy and up 
to relatively large Q'^ » {Ql) (cf. Eq. (4.5)), the calculation of both the inclusive and 
the diffractive cross-sections requires merely the knowledge of the forward amplitude 
{T{b, r))y for the qq pair in the frame where the target carries most of the total energy. 

The dominance of the elastic over the inelastic scattering in diffractive DIS at large 
and relatively small /3 is an essential consequence of fluctuations — in the mean field 
approximation, the opposite situation holds ! [40, 50] — , and it might even seem coun- 
terintuitive at a first sight: At small (3, the projectile contains many gluons, yet Eq. (4.4) 
implies that the relative normalization of these many Fock states remains unchanged after 
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the collision. But in view of the discussion in Sect. 3, this property is in fact natural, as 
we explain now: 

Note first that Eq. (4.4) becomes natural in the saturation regime at r^(Q^) 3> 1, where 
Pdiff and Pel attain their respective unitarity bounds: Pdifr = Pel = 1- In that regime, all 
the Fock states in the projectile are equally absorbed, since they are completely absorbed. 
Thus, for low Q'^ ^ {Q^} at least, the inelastic diffraction is naturally suppressed. Now, 
in the discussion in Sect. 3 we have noticed the tendency of the fluctuations to 'push- 
up the saturation physics', that is, to extend a behaviour which looks natural in the 
saturation regime up to relatively large values of Q^, within a distance ln{Q'^/ (Ql)) ~ o"^ 
above the average saturation momentum which is fixed by the front dispersion. Physically, 
this is so since, within the relevant region at large (the 'diffusive scaling' region), 
the cross-section is controlled by 'black spots' in the target wavefunction, i.e., by rare 
gluon configurations which are already at saturation. Hence, in any scattering event which 
contributes significantly to the diffractive probability (2.3), the projectile with its all Fock 
space components is completely absorbed — the diffraction is purely elastic ! — , and the 
reason where Pdifr is nevertheless smaller than one is just because the probability to find 
such a 'dense spot' is relatively small when Q"^ ^ (Qt)- 

Still from the discussion in Sect. 3, we expect the inelastic probability P^^ff' to be sup- 
pressed with respect to the elastic one P^i by powers of l/cr and/or z/a^, and this will 
be indeed confirmed by the subsequent calculations. But although suppressed in so far as 
the contribution to the cross-section (4.2) — which, we recall, is 'inclusive' with respect 
to the rapidity gap (in the sense of including all the gap values from Ygap to Y) — is 
concerned, the inelastic diffraction represents the leading contribution to the differential 
cross-section per unit rapidity gap at sufficiently small /3 <S 1, and hence cannot be ig- 
nored when computing this particular quantity (even at very high energy). In fact, within 
the present approximations, the contribution of the elastic scattering to the differential 
cross-section d(Tdifr/dln(l//9) is a 5-function peaked aX (3 = 1. This contribution gets 
smeared when going beyond the leading-log approximation with respect to ln(l//3) (see 
Eq. (B.l)), but is in any case locahzed near /3 — 1; hence, for /3 <g 1 we expect most of 
the respective contribution to come from the inelastic diffraction: 

,H,HW) d^/dVi:iC(^.'-;Q)r(-^(6.r;r.w) for „ « 1, 

(4.7) 

where in the r.h.s. Fgap = Y — ln(l//3). 

To summarize, within the high-energy regime specified in Eq. (4.5), the diffractive cross- 
section (4.2) is dominated by its elastic component and can be computed with the qq 
pair alone (cf. Eqs. (4.5)-(4.6)), so like the inclusive cross-section (cf. Eq. (4.3)). But the 
calculation of the inelastic diffraction is still necessary, for two reasons: (i) to explicitly 
check the dominance of the elastic over the inelastic diffraction, and (ii) to compute 
the differential diffractive cross-section (4.7) at small /3 ^ 1. Within the present ap- 
proximations, P^iff^ is given by Eqs. (2.18) and (2.3) which, we recall, are valid so long 
as the saturation effects in the projectile wavefunction remain negligible (this requires 
Q;sln(l//3) <^ ln(l/Q;^), cf. Eq. (2.13)). For any /3 satisfying this condition, Eq. (2.3) can 



41 



be in principle computed by using the dipole probabilities Pn obtained from Monte-Carlo 
simulations of the dipole picture [11] together with the dipole scattering amplitudes con- 
structed in Sect. 3 (or, more generally, obtained by numerically solving the Pomeron loop 
equations [14]). 

But in order to illustrate the main points of physics, it is preferable to use analytic calcu- 
lations, and this is what we shall do in what follows. For these calculations to be tractable, 
we shall restrict our computation of P^^^ to the simplest non-trivial case — that where the 
onium wavefunction at the time of scattering contains only two components: a one-dipole 
state representing the original qq pair, and a two-dipole state corresponding to a qqg con- 
figuration. This is the physical situation in a frame in which the rapidity Y — Yq = ln(l//?) 
of the projectile is relatively large, so its high-energy evolution cannot be neglected, but 
not too large, so that one can restrict oneself to a single step in this evolution: the emission 
of one gluon with longitudinal momentum fraction (3. This restriction on /3 does not entail 
any serious loss of generality in so far as the proof of Eq. (4.5) is concerned — as previ- 
ously explained, the dominance of elastic over inelastic diffraction at high energy follows 
from general arguments, and our subsequent calculation is merely intended to illustrate 
such arguments and to render them more quantitative. On the other hand, this restriction 
would become of problem for phenomenological studies of diffraction at very small values 
of j3; but in that case, it should be straightforward — at least, at conceptual level — to 
extend our analysis by using numerical simulations of the dipole picture within Eq. (2.3). 
Note that the qq and the qqg components have been already included in previous studies 
of DIS in the presence of saturation [36,37,39-43,45-47,49,50,52-54], which were how- 
ever based on the mean field approximation for the target correlations. The comparison 
between such previous studies in the literature and our subsequent results will make it 
easier to emphasize the role of the gluon-number fiuctuations at high energy. 

To conclude this subsection, let us summarize here, for more clarity, the strategy that 
we shall adopt in our analysis: In what follows, we shall evaluate the cross-sections in 
Eqs. (4.3), (4.6) and (4.7) (the latter, at the level of the qqg Fock state) by using the high- 
energy estimates (3.13)-(3.14) for the dipole amplitudes. Note that, whereas the functional 
form of these amplitudes is universal (in the sense explained in the Introduction), the 
coefficients A and which implicitly enter their expressions (via the F-dependencies of 
{ps) and a^) are not: they depend upon the details of the QCD evolution, and presently 
they are not well under control, except in the limit — > (cf. Sect. 3.2). To avoid 
specifying these parameters, we shall perform our subsequent analysis in terms of the 
variables p — (ps) and a. That is, we shall increase the energy by increasing the value of 
a, and we shall measure all dimensionful quantities, like or Q^, in units of (Qs)- This 
strategy is physically meaningful, since the proper way to approach the "high-energy 
limit" is to simultaneously increase Y and Q^, in such a way that the ratio Q^/{Ql) 
remains within the interesting scaling region in Fig. 1. 

To better emphasize the effects of the fluctuations, we shall also estimate the relevant 
convolutions with the following 'mean-field' dipole amplitude 




for r < l/Qs{Y) 
for r > l/Qs(Y), 
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which is representative for cither the Golec-Biernat and Wiisthoff (GBW) 'saturation 
model' [40] (in that case, 7 = 1), or for the solution to the BK equation in the regime 
of geometric scaling (when 7 = 7q 0.63). More generally, for our present purposes, the 
amplitude (4.8) is also representative for the actual situation in QCD at intermediate en- 
ergies, that is, so long as the dispersion of the fronts remains neghgible (of. Eq. (3.20)). To 
study the transition between this regime and the high-energy regime defined by Eq. (3.11), 
where the amplitude takes its "error-function" form in Eq. (3.13), we shall perform nu- 
merical simulations with the dipole amplitude in Eq. (3.17). Indeed, the latter provides 
a smooth interpolation between the mean-field amplitude (4.8) at low energy (cr 0) 
and the error-function amplitude (3.13) at high energy {a — > 00). In most of the explicit 
calculations below, we shall choose 7 = 1, for simplicity. 

To obtain order-of-magnitude estimates in what follows, we shall rely on simplified ver- 
sions of the convolutions appearing in Eqs. (4.1)-(4.2), which in addition to being simpler, 
have also the merit to clearly separate the physical origin of the various contributions. 
The relevant simplifications exploit the properties of the "virtual photon wavefunctions" 
l^/'yy^P, as explained in the Appendix A. The ensuing, simplified, formulae read then 



/ 



f 



2/Q . oc 



dr 



Q2 J 

2/Q 



, (4.9) 



and, respectively. 



2/Q -] °° A 

ldrrnr) + - J ^ f{r 

^ 2/Q 



, (4.10) 



where in writing the r.h.s.'s we have kept the various parametric dependencies, but ig- 
nored all numerical factors. For each of these expressions, the first term within the square 
brackets corresponds to the symmetric configurations, for which v ^ 1/2 (in such config- 
urations, the longitudinal momentum of the incoming photon is "democratically" divided 
among the quark and the antiquark), and the second one, to the aligned jet configura- 
tions, for which v is either close to zero, or close to one (either the quark, or the antiquark, 
carries most of the total longitudinal momentum) . 



4.2 The qq component 



In this subsection, we shall evaluate the DIS cross-sections which can be computed from 
the scattering of the qq alone, namely, the inclusive cross-section (4.3) and the elastic 
component of the diffractive cross-section, Eq. (4.6). As anticipated in the previous dis- 
cussion, Eq. (4.6) is in fact the dominant contribution to the diffractive cross-section (4.2) 
at high energy. Accordingly, the cross-sections that we shall compute in this subsection 
represent our main results in this paper. 

To render the subsequent discussion more intuitive, it is instructive to see the plots of the 
integrands in Eqs. (4.3) and (4.6), as numerically computed with the two limiting expres- 
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Fig. 12. The qq contributions to the integrands in Eqs. (4.3) and (4.6) as computed with two 
different expressions for the dipole ampUtude — the mean-field ampUtude (4.8) with 7 = 1 
and the 'fluctuation' ampUtude (3.13) — are plotted as functions of rQg for two values of : 
Q2 = o.ig2 ^left) and Q"^ = lOOQ^ (right). 



sions for the dipole amplitude: the high-energy approximation (3.13) (the 'fluctuation' 
piece in the plots in Fig. 12) and the mean-field model of Eq. (4.8) (with 7 = 1 for defi- 
niteness). As manifest on Fig. 12, and it will be analytically demonstrated in what follows, 
the strength of the integration is shifted towards smaller dipoles sizes after including the 
effects of fiuctuations. This phenomenon becomes truly spectacular at high Q^, where the 
integrand corresponding to 'fluctuations' appears to be strongly peaked at r ~ 2/Q (for 
both inclusive and diffractive cross-sections), while that for the mean field approximation 
is rather smoothly distributed at all sizes r ^ 2/(5 (with only mild maxima at r ~ 2/(5 
for (Jtot and, respectively, at r ~ l/(5s for (Jdifr)- 



4.2.1 The mean-field approximation 

Let us first present the corresponding mean-field estimates, which are well-known (see, 
e.g., Refs. [39, 40, 54]), but are interesting here as a term of comparison with the forthcom- 
ing results which will include the effects of fluctuations. These estimates will be obtained 
here by using the dipole amplitude (4.8) together with the approximations in Eqs. (4.9)- 
(4.10). To avoid a proliferation of physical situations, we shall focus on the case ^ (5^, 
where the contrast with the high energy behaviour turns out to be most important. For 
that case, we shall distinguish between the total and the diffractive cross-sections, and 
also between transverse and longitudinal polarizations for 7*. To simplify writing, we shall 
omit the overall factor ctem-ZVc e^, which is common to all the cross-sections. 

i) Total cross-sections 

Choosing 7 = 1 in Eq. (4.8), the transverse contribution to Eq. (4.3) is found as 




(4.11) 
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where the dominant term, as isolated in the r.h.s., comes from the ahgned-jet configura- 
tions with dipole sizes r within the range 2/Q < r < 1/Qs- For the longitudinal sector, 
one similarly obtains 




(4.12) 



where the leading-order term comes from dipole sizes r ~ l/<5, and it receives equally 
important contributions from both the symmetric and the aligned~jet configurations. 

The above results show that, at sufficiently large Q^, the transverse sector dominates over 
the longitudinal one, by a logarithm \n{Q'^/Ql). However, this dominance disappears at 
lower values of Q^, within the geometric scaling window in Fig. 1. Indeed, in that regime 
one must rather use 7 = 70 ~ 0.63, and then one finds that both (7^(. and cr^t are of the 
same order, namely ~ {Ql/Q'^V^^ ^^^Y dominated by dipole sizes r ~ l/Q. 

After similarly treating the case Q^, one finds the following limiting behaviours for 

the total 7*/i cross-section in this mean-field scenario : 

- { 4.13 

[Ql/Q'HQVQD forg2»g2. 



ii) Diffractive cross-sections 

The diffractive cross-section (4.6) involves the dipole amplitude squared, and this leads 
to important differences with respect to the inclusive cross-section, at least in this mean- 
field scenario. The contribution of the small dipoles is now strongly suppressed, because 
~ (r'^Ql)'^^ for r -C l/Qs- Accordingly, the dominant contribution to cxdifr comes from 
relatively large dipole sizes, of the order of the saturation length: r ~ l/Q^ [34, 35, 40, 68]. 

To see this, let us focus on the contribution of the aligned-jet configurations within the 
transverse sector (this turns out to be the dominant piece). The relevant integral reads 




(4.14) 



where for more clarity we have split the integration domain into two: 2/Q < r < 1/Qs 
and r > 1/Qs- It can be easily checked that, so long as 7 > 1/2, both domains contribute 
on equal footing to the leading order result ~ Q^/Q^, which is generated by dipole sizes 
r ~ 1/Qs, as anticipated. Indeed, the integral over the first domain is dominated by large 
values of r, of the order of the upper cutoff 1/Qs, while that over the second domain is 
saturated by small values of r, close to the respective lower cutoff 1/Qs- By comparison, 
the other contributions to diffraction — that of the symmetric transverse configurations, 
and the whole longitudinal contribution — are of 'higher-twist order', that is, they are 
suppressed by powers oi Ql/Q'^ <^ 1 with respect to the leading contribution (4.14). 

We thus arrive at the (by now) standard picture of DIS diffraction in the presence of 
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saturation [34,35,40], in which adm is a leading-twist quantity at high-Q^, so Uke the 
inclusive cross-section (4.11). Accordingly, the ratio between the diffractive and the in- 
clusive cross-sections is only slowly varying with and with the total energy [40, 47] 



" (da.„,/d'(,) \n(QyQi(Y)) ' ^ ' "> 

(this behaviour is illustrated by the 'mean-field' curve in Fig. 17), in rough agreement 
with the pattern observed in the corresponding HERA data at small x. This agreement 
has represented one of the main successes of the saturation models applied to the phe- 
nomenology at HERA [40,47]. 

Notice the qualitative difi^erence between the ways that the inclusive and diffractive cross- 
sections get constructed — (Ttot receives leading-twist contributions from all dipole sizes, 
with a slight preference though for r ~ l/Q, whereas adifr is dominated by dipoles with 
r ~ l/Qs — , which nevertheless leads to quantitatively sirailsn results (cf. Eq. (4.15)). For 
later reference, let us display here the qq contribution to the diffractive cross-section in 
the mean-field approximation for both small and large Q^: 



da'" 



diflf 



HQl/Q') for g2 « Q2 

(4.16) 

Ql/Q' for » . 



The above results depend crucially upon the property of the mean-field amplitude T(r) 
to decrease very fast when decreasing r below l/Qg : T ^ (^^Qs)^ with 7 > 1/2. From the 
discussion in Sect. 3, we know that this behavior is eventually washed out by fluctuations 
when increasing energy. As we shall shortly see, this change of behaviour has dramatic 
consequences for the high-energy hmit of both inclusive and diffractive DIS cross-sections. 



4.2.2 The high-energy behaviour: 'black spots' and diffusive scaling 

We now turn to the most interesting case for us here, namely the high-energy regime of 
Eq. (3.11). The kinematical region where the effects of the fluctuations are most visible, 
and that we shall concentrate on in what follows, is the high-Q^ range deflned by 

a(F)>l/7o and a{Y) < Z = In ^ < a^{Y) . (4.17) 

{Qi) 

For much smaller Q^, such that <^ Ql{Y), the cross-sections are dominated by 'sat- 
urated' dipoles with sizes r ^ for which the mean fleld approximation works quite 
well. For much larger Q^, such that Z ^ cr^(^), one enters the standard perturbative 
regime where the amplitudes show color transparency and evolve according to the DGLAP 
equation. 

We shall soon verify that, within the interesting range (4.17), the convolutions yielding 
the DIS cross-sections are dominated by dipole sizes r ~ 1 /Q (at variance with the mean 
fleld case, especially in the case of diffraction). It is therefore appropriate to evaluate 
these convolutions with the expression (3.13) for the dipole amplitude, which is valid for 
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z -C By using that, wc shall be able to compute exactly the dominant behaviour 

of the cross-sections at high energy. But before describing the exact calculations, we 
shall present some parametric estimates which have the merit to simply demonstrate the 
physical points that we would like to emphasize here. 



i) Total cross-sections 

We start with the transverse sector and, for more clarity, we separate the respective 
contributions of the symmetric and the aligned-jet configurations: 

^ - SL + , (4.18) 
where (with z = ln{l/r^{Ql)) ) : 

T r dr f°° , „ „ f z 



/ * {nr))y ~ /~dz Eric (i) , (4.19) 



^L-^J'^ mr))y ~ £ d. e' Erfc (i) . (4.20) 

2/Q 

Consider the large-Q^ case, such that a <^ Z <^ a^. It is then easy to verify that: (a) 
the integrals giving 5*^^. and are dominated by their respective endpoints at z — Z 
(i.e., by r ~ 1/Q), (b) the dominant behaviour at large Z is given by the Gaussian 
exp(— ^^/o"^), and (c) after also computing the prefactors in front of this Gaussian, the 
dominant contribution is found to come from the symmetric configurations. 

To demonstrate this, we shall first consider the symmetric integral: 
SLt ^ a I dx Erfc(x) a i — e a — e 

Jz/a Jz/a X Jzya'' t 

where x — z/a, t — x^, and we have used the asymptotic behaviour of Erfc(x), cf. 
Eq. (3.18) (recall that Z/a ^ 1). It is now clear that the integral over t is dominated by 
its lower end at t — Z'^/a^, i.e., hy z — Z, as anticipated. 

The aligned-jet contribution can be similarly treated: because of the exponential factor e^ 
within the integrand, the integral in Eq. (4.20) is dominated by its upper cutoff at z — Z, 
which yields 




(4.22) 



As anticipated, this contribution is suppressed (by a factor Z/cr^ -C 1) with respect to 
the contribution (4.21) of the symmetric configurations. 
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Similar conclusions hold for the longitudinal sector : the dominant contribution is again 
generated by dipole sizes r ~ l/Q and it is of the same order as A^t? Eq. (4.22), hence it is 
parametrically suppressed. We thus conclude that the inclusive cross-section in the high- 
energy limit and for relatively large is dominated by symmetric dipole configurations 
with typical sizes r ~ l/Q within the transverse sector. 

It is in fact quite straightforward to go beyond the parametric estimate in Eq. (4.21) and 
compute the dominant behaviour exactly. We present here only the final results, since a 
similar calculation will be explained in more detail below, for the case of diffraction. Specif- 
ically, after also reintroducing the overall electromagnetic factor F = (iVcaem/27r^) J^fc'j, 
one finds that the longitudinal part of the total cross-section is given by 

d^-.^fErfcfiV (4.23) 



d% 3 \a 
while the respective transverse part reads 



"^""^ (4.24) 



d% 3 '\a 
Here we have found convenient to define the function 



$i(x) = /dt'Erfc(i;) = ^ exp(-a;^) - a;Erfc(a;) (4-25) 

J -v/tT 

X 

which has the following behavior in the various limits 

;xp(— 
2^x^ 



2k| + ^?fc$ for.«-l 



1 

exp(— x^) 
20Fx2 



for X = (4.26) 
for a; 3> 1. 



Eqs. (4.23) and (4.24) hold for Z <^a^ and are correct up to corrections of relative order 
Zja^ and/or \ ja. They confirm the previous estimates hke (4.21) and imply the following, 
final, result for the total DIS cross-section in the high-energy regime: 

Note that the quantity (d(Ttot/d^6)/cr shows diffusive scaling, i.e., it depends upon the 
kinematical variables and Y only via the dimensionless variable Z/a = ln((5^/ {Ql))/a. 

Eq. (4.27) also allows us to study the behaviour of the cross-section at relatively low 
(below and around (Ql)). Deeply inside the saturation region, one finds 

^"^^"^ ln((Q^)/g^) »a, (4.28) 

which is in agreement, as expected, with the mean-field behaviour in Eq. (4.13). However, 
the difference with respect to the mean field situation occurs already in the vicinity of 
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Fig. 13. The total cross section as a function of Q"^/ {Ql) and for various values of the front 
dispersion a. 
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Fig. 14. The total cross section (divided by the front dispersion) as a function of the "dif- 
fusive" scaling variable Z = ln{Q'^ / {Qg)) / a . 



the (average) saturation line Z = : Within a wide interval | ln{Q'^ / {Ql))\ ^ a around 
this line, the cross-section (4.27) is rather large — since proportional to cr — and Q^- 
independent: 



a for |ln(gV(g^))| « ^. 



(4.29) 
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The cross-section starts to fall only when Q becomes so large that ln(Q7(Q2)) > a, 
and then it has a rapid fall, of the Gaussian type (cf. Eq. (4.21)). 

This behaviour is illustrated in Figs. 13 and 14, which are the analog at the level of the 
DIS inclusive cross-section of Figs. 10 and 11 for the dipole amplitude. Namely, these 
figures represent results obtained via the numerical evaluation of the total cross-section 
(4.3) with the dipole amplitude in Eq. (3.17) with 7 = 1. Fig. 13 shows the increase 
in the deviation from the mean-field behaviour (and thus from geometric scaling) when 
increasing a, whereas Fig. 14 demonstrates the emergence of the diffusive scaling and the 
convergence of the cross-section towards its high-energy asymptotic in Eq. (4.27). 

ii) Diffractive cross-sections 

Moving to the diffractive sector, one can immediately notice an important difference with 
respect to the mean-field scenario of Sect. 4.2.1: Whereas in that case, the replacement of 
the dipole amplitude by its square has strongly suppressed the small-dipole configurations 
and thus shifted the strength of the integration from r~ 1/Q to r ^ l/Qs, in the present, 
high-energy, case, there is no similar suppression, because both the Gaussian exp(— ^^/cr^) 
and its square exp{—2z'^/a'^) decay rather slowly (as compared to the exponential e~^) 
within the interesting range at a <^ z <^ a^. Accordingly, the aligned-jet contribution to 
high-energy qq diffraction, namely (compare to Eq. (4.14)) 

Q./^TO>J~e-/_^d.e'Erfc^(i). (4.30) 

2/Q 

is dominated by its lower end at r ~ l/Q (that is, by 2; = Z), so like the corresponding 
contribution (4.20) to the inclusive cross-section. 

This points towards an important physical difference between the mean-field (or inter- 
mediate energy) behaviour and the behaviour at high energy, which deserves a more 
qualitative explanation. To that aim, notice that the measure dr/r^ in the aligned-jet in- 
tegral favors the small dipoles, so the integral would be naturally dominated by its lower 
end at r ~ l/Q ii there was not for the strong suppression of the dipole sizes r <^ 1/Qs 
introduced by the dipole amplitude squared T^. Such a strong suppression occurs, as we 
have seen, for the mean-field amplitude T, and also for the typical events in the statisti- 
cal ensemble in QCD at high energy and large ^ (Ql) : indeed, a typical event has 
Qs ~ (Qs) and thus yields a small contribution, of order (Q^) /Q^ = e^^, to the inte- 
gral. However, due to the front dispersion in the presence of fluctuations, the statistical 
ensemble contains also fronts which are at saturation at the minimal size l/Q, and each 
such a front yields a relatively large contribution, of 0{1). Although such fluctuations are 
relatively rare, their contribution weighted by the respective probability ~ exp(— Z^/o"^) 
is still larger than that of the typical fronts, which behaves like c^^. In other terms, 
the convolution peaks up those rare gluon conflgurations, or 'black spots', in the target 
wavef unction which are at saturation at the scale set by the virtual photon. 

Therefore, at high energy — and in contrast to what happens in the mean-fleld approxi- 
mation, or at intermediate energies — the inclusive and diffractive cross-sections receive 
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their dominant contributions from the same physical configurations, namely the small 
dipoles with size r ~ In fact, a qualitative analysis similar to that in Eqs. (4.18)- 
(4.22) reveals that the analogy between inclusive and diffractive processes at high energy 
is even stronger: In the large-Q^ regime defined by Eq. (4.17), both types of cross-sections 
are dominated by symmetric dipole configurations within the transverse sector. In par- 
ticular, the dominant behaviour of cr^fg can be estimated as in Eq. (4.21), and reads 



^ si, ~ ^3 exp (-^ j for a«Z«a\ (4.31) 

which implies the following estimate for the (diffractive/inclusive) ratio R : 

i? ~ exp (^-^j for (7 < Z < (7^ (4.32) 

This is a scaling function, which decreases quite rapidly with Q"^ at fixed Y (unlike the 
respective mean field estimate in Eq. (4.15)), but one should notice that this decrease 
becomes significant only at very large Q^, such that Z ^ a. 

To make these estimates more precise, let us compute the dominant behaviour of a'^fg 
in the regime ajQ ^ 1 and Z a^. This is obtained by evaluating the convolution in 
Eq. (4.6) with the dipole amplitude (T(r)) in Eq. (3.13). Consider first the longitudinal 
cross section. Using Eq. (A. 2) and defining u — rQ and r — Q/ (Qs) we have 

d^difr _ ttF /• 7^^_ ..^ / ^ /— -\ ^ . /ln(rVM2)' 



^ J dt- J du Av^{l - u Kl [u^v{l-v)^ Erfc^ ( ) " (^-^S) 



When (7 is large we can set rt = 1 in the argument of the error function; indeed, as 
previously explained, the convolution is dominated by r ~ l/Q- Then the integrals over 
u and V are easily performed and we arrive at (note that Inr^ = ln((5^/((5s)) — Z) 

Now let us look at the transverse cross section. It reads 

^ ^ ^yd^Jdt.^(i-^)[^2+(i-^)>K2(z,y;;(r^) Erfc^ Ml^Y 

^ ^ 

(4.35) 

The behavior of the last factor in the small-w region is crucial since it cancels the loga- 
rithmic singularity of the remaining part of the integrand at -u = 0. Hence, we cannot set 
M = 1 as we did for the longitudinal case. Here it is convenient to define the function 

oo I — 

$2(2;) = / dvETic^{v) = exp(-x2) Erfc(a;) - xErfc^(a;) - J- Erfc(\/2a;), (4.36) 

J \/tT V TT 

X 

which has the following behavior in the various limits 
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for X <^ —1 



for X — 



for X ^ 1. 



(4.37) 



Using Eq. (4.36) we can write the last factor in Eq. (4.35) as 



Erfc^ 



2"d^^^ 



(4.38) 



and then the li-integration in Eq. (4.35) can be performed by parts. The boundary term 
vanishes and therefore the virtual photon transverse cross section becomes 



1 T 1-1 -'- ^ 



v)f^ [v^ + (1 
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X Ko (u^Jv{l - v)j Ki (^u^Jv{l - v) 



(4.39) 



Now we can safely set w = 1 in the argument of the $2 function since the remaining part 
of the integrand is well-defined for any value of u. The integration over u and v becomes 
straightforward and we obtain 



a $2 - 
6 Vcr 



(4.40) 



Comparing Eqs. (4.34) and (4.40) we see that the transverse cross section dominates over 
the longitudinal one within the whole region of interest (namely, — 00 < Z <^ cr^ and 
(T7o ^ 1), so we finally arrive at 



dcrX 



diflf 



dal 



(4.41) 



up to corrections of relative order Z/a'^ and/or 1/cr. As anticipated by the notation in 
Eq. (4.41), this is in fact the complete dominant contribution to the diffractive cross- 
section in the high-energy regime of interest (see the discussion in Sect. 4.3). This quantity 
has the same qualitative behaviour as the inclusive cross-section in Eq. (4.27), and this 
behaviour is illustrated in Figs. 15 and 16 (to be read by analogy with Figs. 13 and 14). 
In particular, the rescaled quantity {da2^g/d'^b)/a exhibits diffusive scaling. 

Finally, the ratio R of the diffractive to the total cross section is easily obtained from 
Eqs. (4.41) and (4.27). It reads 



l$2(ln(QV(g.)^)/a) 
2 $1 (ln(gV {Q,f)/a) ' 



(4.42) 



For small values of Q such that ln{{Ql)/Q'^) ^ a the ratio approaches 1, it is constant and 
equal to 1 — \/2/2 = 0.293 in the wide interval | ln((5^/ {Ql))\ <C c around the saturation 
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Fig. 15. The contribution of the qq component of the onium wavefunction to the diffractive cross 
section as a function of Q'^ /{Ql) and for various values of the front dispersion. 




line, while it falls very fast with increasing momentum when \n{Q'^ / (Q'j.)) ^ a with the 
precise form easily inferred from Eqs. (4.37) and (4.26). This behaviour is graphically 
illustrated in Fig. 17, and compared to the corresponding mean-field behaviour. 
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Fig. 17. Tie ratio of the diffractive (contribution of the qq component) to the total cross section 
as a function ofQ^/{Ql). 



4-3 The qqg component 



In this subsection, we shall analyze the first non-trivial contribution to the inelastic 
diffraction, which arises when the wavefunction of the projectile at the time of scattering 
is a superposition of two Fock space components: the original qq pair (one dipole) and 
its qqg excitation (i.e., two dipoles at large Nc). This configuration is physically relevant 
for relatively large diffractive mass Mjr ^ Q^, oi small values of /3 -C 1, but such that 
Q;sln(l//9) <S 1, in order for the effects of the higher Fock states in the projectile to be 
negligible. As we shall see, the two-component configuration generates contributions of 
order ln(l//3) to both the elastic and the inelastic pieces of the diffractive cross-section, 
but at high energy, the inelastic piece is parametrically suppressed with respect to the 
elastic one in the whole kinematical range of interest (including the large-Q^ domain 
defined by Eq. (4.17), which extends with Y). This means that, if one starts at y = Fgap, 
where the diffraction is purely elastic (cf. Eq. (2.10)), and one increases Y for a fixed value 
of Ygap (i.e., one increases the rapidity of the projectile), then the elastic piece of diffraction 
will rise faster than the inelastic one, and thus will remain the dominant contribution to 
(Tdiff at any Y > Yg^p (i.e., for arbitrarily small /?). This situation occurs for sufficiently 
high values of Yg^p, such that cr( Ygap) 3> 1, and is in sharp contrast with the corresponding 
mean-field scenario [43, 50], where for Q'^ » the inelastic diffraction rises faster than 
the elastic one with 1//3 and thus rapidly dominates over the latter. 

One should also remind here that the diffractive cross-section that we have in mind, cf. 
Eqs. (2.1) and (4.2), is integrated over all the values of the rapidity gap between Ygap and 
Y. But, of course, the (iiij^^eren^m/ cross-section dcrdifr/dln(l//3) at sufficiently small /? ^ 1 



Strictly speaking, our argument applies so long as ln(l//?) <C ln(l/a^), i.e., for a projectile 
which is non-saturated, since this is the validity limit of Eq. (2.3). 
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is always dominated by inelastic scattering (cf. Eq. (4.2)), because the corresponding 
elastic contribution is peaked near (3 = 1. In what follows, we shall compute the lowest 
order contribution to this differential cross-section, as given by the qqg state. 



4-3.1 Diffractive DIS with one or two dipoles 

Within the present approximations, the scattering probabilities for an onium made with 
two dipoles are obtained by taking iV = 2 in the general formulas in Sect. 2 and using the 
dipole probabilities Pi and P2 given by Eq. (2.14) with dY = Y — Yq = ln(l//5). 

We are primarily interested in diffraction, but in order to separate its elastic and inelastic 
components, we also need the forward amplitude A{x, y; Y), so let us start with this latter. 
According to Eq. (2.9), the forward amplitude can be expressed either as the scattering 
of a single qq dipole in the frame where the target has rapidity Y : 

A{x,y;Y) = {T{x,y))y, (4.43) 

or as the scattering of an onium composed of two dipoles, in the frame where the target 
has rapidity Yq (with the simplified notation Txy = T{x,y)) : 

Aix,y-Y) = {T{x,y))Y„ + ln{l/ f3) A A{x,y;Yo) , 
AA{x, y; Yo) = ^ J M{x, y, z) (T,, + - T^^T^y - T,^)^^ . (4.44) 

z 

The quantity describes the changes in A brought in by the one-step evolution, namely, 
the scattering of the two-dipole state (the first three terms within the brackets) together 
with the 'virtual' correction to the scattering of the one-dipole state (the fourth term 
there). By comparing Eqs. (4.43) and (4.44) and taking the limit dY = ln(l//3) — > 0, one 
can immediately deduce an evolution equation for {Txy)Y '■ 

= ^jM(x,y, z) {T^^ + T^y - T^y - Tx.T,y)y , (4.45) 

z 

which is recognized as the first equation in both the Balitsky hierarchy [5] and in the 
hierarchy of equations with Pomeron loops [14]. (The differences between these two hi- 
erarchies, which express the effects of gluon number fiuctuations, start with the second 
equation, as obeyed by the two-dipole amplitude {T{xi, yi)T{x2, y2))y ) Within the mean 
field approximation (2.32), Eq. (4.45) reduces to the BK equation (2.39). 

Consider similarly the elastic probability, Eq. (2.7). This too can be evaluated in either 
frame (the 'F-frame' and the 'lo^frame'), with the following results : 

P.i{x,y;Y) = \{T{x,y))Y\' ^ \{T{x,y))yf + (5) AP,,{x,y;Y^) 
AP,,(x, y- Fo) = 2 Re{{T{x, y))Y, AA{x, y; Yo)} , (4.46) 

where the expression in the lo^frame holds to the order of interest in cxg ln(l//3). Clearly, 
APei represents the qqg contribution to the elastic probability. 
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For the diffractive probability (2.3), on the other hand, the choice of the Yo^frame is 
mandatory, since Yq plays also the role of the physical rapidity gap. We thus find: 



P^s{x,y;Y,Yo) = P!Ux,y;Yo) + \n{l/ (3) APa^{x,y;Yo), 
where 



P!Ux,y;Y,) = \{T{x,y)) 



Yo 



Peiix,y;Yo) 



(4.47) 



(4.48) 



is the contribution of the elementary qq pair, and is the same as the elastic probability 
corresponding to a total rapidity Yq (cf. Eq. (4.46)), while 

APMx,y;Yo) ^^J M{x,y,z){\{l - S^,S,y)Yf - |(1 - S^y)yf} 

z 



represents the additional contribution due to the qqg state. This contains both elastic 
and inelastic components, and in order to separate the latter we use P^^g^ — PdiS — Pel, 
together with Eqs. (4.46) and (4.48). This yields: 



P^ix, y- F, Fo) = ln(l//9) AP.U^, y, Yq) 
An^,,{x,y;Yo) = AP^^{x,y-Yo) - AP,,{x,y;Yo) . 

Prom the expressions (4.44), (4.46) and (4.49), one immediately finds 



APinei(£c, y; Yq) 



as 
2^ 



M{x, y, z) (T,, + T,y - T^y - T^^Z 



(4.50) 



(4.51) 



Note that the inelastic diffraction starts at order Q;sln(l//9), i.e., it requires at least one 
soft gluon in the wavef unction of the projectile. 

From Eqs. (4.7) and (4.50), one can deduce our lowest-order approximation for the dif- 
ferential diffractive cross-section, valid for /3 <^ 1 but such that agl'o.il/ j3) -C 1 : 



^ff—-^['dz[d'r ^ \ra{z,r;Q)\' APUb,r;YQ). 
ln(l//3) Jo J 



(4.52) 



Note that, for the simple onium configuration at hand, the inelastic diffraction is, by 
construction, the probability to have two dipoles in the final state. Hence, the quantity 
(4.52) represents also the cross-section for gluon production in diffractive DIS [45,52], at 
the present level of accuracy. 



4-3.2 The inelastic component of diffraction 

The inelastic diffraction probability associated with the qqg state is given by Eq. (4.51), 
which after the trivial change of variable z = a; — s is rewritten as (with r — x — y) : 
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APinei(r) = ^ / /"_ (T(5) + T(r - - r(r) - T(s)T(r - s)f . (4.53) 

In line with the previous approximations, we have assumed the high-energy amphtudes to 
be real and ignored their impact parameter dependence. It is also implicit in our notations 
that the target expectation values are evaluated at a rapidity Yq, and therefore they involve 
the (average) saturation momentum {Ql) (Yq) and the front dispersion cr'^iYo) = Dfj-agYo. 
We shall evaluate this expression in both the mean-field approximation and the general 
case with fiuctuations, to emphasize the difference between these two situations. 



i) The mean-field approximation 

In this approximation, the two-dipole amplitude is assumed to factorize : (T(s)T(r — s)) ^ 
T{s) T{r — s), with T(r) the mean-field amplitude of Eq. (4.8). The most interesting situ- 
ation is when the external dipole is small, r -C ^/Qs, and for that case one can distinguish 
between three different physical regions for the internal dipoles s and r — s : (i) one of 
the two dipoles, say s, is much smaller than the external one: s <S r; (ii) both internal 
dipoles are much larger than the external one, but they are still small as compared to the 
saturation length: r ^ s ^ |r — s| ^ ^/Qs] (iii) both internal dipoles are at satura- 
tion: s ~ |r — s| ^ ^/Qs- Corresponding to these three cases, the integral in Eq. (4.53) is 
decomposed into three pieces, which are estimated as (up to irrelevant numerical factors) 



APi„ei(r) ~ as 



r ^/Qs oo 

r i/Q^ 



ds 



asr'Qi (4.54) 



The dominant contribution, as isolated in the r.h.s., comes from relatively large internal 
dipoles with sizes s ~ l/Qs ^ r (i.e., from domains (ii) and (iii) alluded to above). 
Note the similarity between this calculation and that of the qq contribution to the diffrac- 
tive cross-section in the MFA (cf. Eq. (4.14)) : the external dipole size r plays here the 
same role as the resolution scale 1/Q in Eq. (4.14), and the dipole kernel has limiting 
behaviours similar to those of the transverse virtual photon wavefunction (cf. Eq. (4.9)). 

The case of a relatively large dipole r ^ 1/Qs is technically more involved, but the 
analysis in Ref. [50] shows that, in that case, APinei(r) vanishes very fast — at least as 
fast as S'^{r) = (1 — T{r)Y — when increasing r. 

To summarize, in the mean field approximation, APinei(r) increases with r like at 
r^Ql ^ 1' it decreases very fast at r'^Ql S> 1, and it develops a maximum at r ~ 1/Qs 
[50]. This behaviour is manifest on the "mean field" curve in Fig. 18. 



ii) The high-energy behaviour with fluctuations 

We shall now extract the dominant behaviour of the integral in Eq. (4.53) in the high- 
energy regime at a 3> I/70 and z = \Yi{l/r'^{Q^g)) <^ a^. We anticipate that, in this 
regime, the dominant contribution comes from internal dipoles with sizes comparable to 
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that of the external one, s ~ |r — s| ~ r, hence one can use the high-energy estimates 
(3.13)-(3.14) for all the amplitudes which appear in Eq. (4.53). 

Eq. (3.14) implies (T(s)T(r — s)) = (T(r<)), where r< = min(s, |r — s|); thus, the 
quadratic and one of the linear terms in Eq. (4.53) cancel each other, and the integral 
reduces to 



APinel — 2 ■ 



2 



2% J S^Cy. _ 



l2 



{T{s)) - {T{r)) , (4.55) 



a 



where the integration is now restricted to the half plane \s\ > \r — s\. It is convenient to 
change the integration variable to u — s/r. Then the above equation becomes 



1/2 



where the upper limit is 0o — arccos(l/2M). Since the scattering amplitudes in Eq. (4.56) 
do not depend on the angle 0, one can perform the angular integration to obtain 



1/2 ^ ^' ' 

One can check that the above integrand vanishes as for large u. Thus, the dominant 
contribution comes from u — 0{1), that is, from internal dipoles sizes s such that s ~ r, 
as anticipated. By using this property together with the specific form of the amplitude 
given in Eq. (3.13), we can expand 

(T(ur)) - (T(r)) . -— .^(.^(j^^^)^^) = exp [-^^ ) . (4.58) 

where we have kept only the first term in the Taylor expansion, since the other terms 
are suppressed by higher powers of ln{u^)/a and thus are truly negligible at high energy. 
After the expansion (4.58), the result of the integration is a pure number, so we arrive at 

APinei(r) ~ ^ ^ exp for z « (J^ (4.59) 

with J numerically computed as J = 15.1. 

The above result is valid for both small and large dipoles and shows that APinei(r) is 
an even function of 2; = \n{l / r'^ (Ql)) , which rapidly vanishes when \z\ ^ cr and has a 
maximum at r'^{Q^) = 1. This behaviour is compared to the corresponding mean field 
behaviour in Fig. 18 which shows that, when increasing the energy, the maximum around 
~ l/(Qs) becomes flatter and flatter. The emergence of the diffusive scaling with 
increasing a and the convergence towards the asymptotic form (4.59) are illustrated in 
Fig. 19. (The various curves in Figs. 18 and 19 are obtained by numerically evaluating 
the r.h.s. of Eq. (4.53) with the dipole amplitude in Eq. (3.17) with 7 = 1, together with 
a corresponding expression for the two-dipole amplitude (T(s)T(r — s)).) 



(l+tt) / 2m-1 \ 
l-u| V 2m+1 / 
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Fig. 18. The contribution of the qqg component of the onium wavefunction to the inelastic 
diffraction probabiUty as a function of r^{Q^) and for various values of the front dispersion. 
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Fig. 19. The same as in Fig. 18 (multiplied by ~ cr^j as a function of z/a. 



By inserting the above result (4.59) into the r.h.s. of Eq. (4.52) and performing the con- 
volutions with the photon wavef unctions, we can compute the lowest-order contribution 
to the differential diffractive cross-section at high energy and /3 <S 1. (As previously 
explained, this quantity is the same as the the cross-section for gluon production in 
diffractive DIS.) Via manipulations similar to those encountered before (cf. Sect. 4.2.2), 
one finds 
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da^iff as 2JF 1 f 2Z 



— exp , (4.60) 



for the longitudinal piece and, respectively, 

^-"-^ ^^^lErfcf^^ (4 61) 

d26dln(l//?) - 27r 6 a [ a )' ^ ' 

for the transverse one. As usual, F = [N centra/ 27^"^) J2f^}: and the above results are valid 
in the high energy regime at a^o » 1 and — oo < Z <^ a^, up to corrections of 0{Z/a^). 

Comparing Eqs. (4.60) and (4.61) we see that the transverse cross section dominates the 
longitudinal one everywhere in the region of interest. Therefore, 



d26dln(l//3) 2% Q a 

Recall that the average saturation scale (Ql) and the front dispersion cr^ which enter 
this equation is to be evaluated at a rapidity Yq = Vgap- The cross-section (4.62) is a 
monotonically decreasing function of (for fixed Y, that is, fixed a and {Ql))- Note 
that after multiplication by a, this cross-section exhibits diffusive scaling. 

It might be interesting to compare this high-energy estimate to the corresponding result 
at intermediate energies, as given by the mean field approximation (say, by the solution 
to the BK equation). In the high-Q^ regime at ^ Ql, the latter can be estimated by 
using APinei ~ CKs^^Qs, cf. Eq. (4.54), which yields [50] 

"^^^'^^ '^sF%\n% for g^»g^ (4.63) 



d26dln(l/^) * Q 

This mean-field result receives contributions from dipole sizes r uniformly distributed 
within the range 1/Q < r < 1/Qs- By contrast, the high-energy estimate (4.62) is domi- 
nated by r ~ l/Q- 



4-3.3 The elastic diffraction 

In this subsection, we shall evaluate the increase APci in the probability for elastic scatter- 
ing due to the emission of one soft gluon (cf. Eq. (4.46)). Our purpose is to demonstrate 
that, in the fluctuation-dominated regime at high energy and relatively large Q^, cf. 
Eq. (4.17), this increase is parametrically larger than the corresponding evolution APinei 
in the probability for inelastic diffraction, that we have already computed. 

Let us start by briefly discussing the case of the mean field approximation, where the 
opposite situation holds: for Q'^ ^ Ql (or relatively small dipole sizes), the inelastic 
piece dominates over the elastic one. Indeed, by using APei(^) = 2 {T{r)) AA{r) together 
with A^ = d {T)/dY, cf. Eqs. (4.44)-(4.45), and the expression (4.8) for the mean-field 
amplitude, one immediately finds: APei ~ ag (r^Qs)^^ for r <^ l/Qs- This is a 'higher- 
twist' effect, much smaller than the corresponding inelastic contribution in Eq. (4.54). 
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But at high energy, APci turns out to be the dominant contribution, down to very small 
dipole sizes, namely, so long as z <^ a^. Moreover, its calculation turns out to be quite 
subtle, as we explain now. On one hand, it is straightforward to use = d {T{r)) / dY 
together with Eq. (3.26), to deduce : 



which is indeed parametrically larger (for any z) than the corresponding inelastic piece, 
Eq. (4.59). For instance, for small dipole sizes, or large z ^ a, this elastic component 
behaves like 



and dominates over the inelastic piece by a large factor Xa^/z ^ 1, as anticipated. 

On the other hand, the above estimate for APci has the drawback to involve the average 
front velocity A, which is generally unknown. It is therefore tempting to try and compute 
AA{r) directly from its integral representation (4.44) together with Eqs. (3.13)-(3.14) 
for the dipole amplitudes. A priori, this calculation is independent of A and should even 
enable us to determine A, by comparison with the expected result in Eq. (3.26). But this 
expectation appears to be naive: the coefficient of the would-be dominant term at large Y 
— that term precisely which should be identified as A/o" — turns out to be zero. (This can 
be checked via manipulations similar to those in Eqs. (4.56)-(4.58): without taking the 
square of the amplitudes inside the integrand, so like in Eq. (4.57), the integral vanishes 
exactly after the expansion (4.58).) This seems to imply the unphysical result A = 0, but 
in reality it only means that the approximations (3.13)-(3.14) — although the correct 
leading-order estimates for the average amplitudes at high energy — cannot be used to 
also evaluate the right hand sides of the evolution equations. The latter are sensitive to 
the subleading terms, which are not under control in the present approximations. 

To summarize, in the high-energy regime of interest, the diffractive probability associated 
with the qqg component is controlled by its elastic piece, APjifr ~ APgi, and has the 
dominant behaviour displayed in Eq. (4.64). Although not as symmetric as the inelastic 
contribution (4.59), the function (4.64) has still the properties to vanish for both very 
small and very large dipole sizes, and to develop a rather flat maximum at r ~ 1/ {Qs)- 
This behaviour is illustrated in Fig. 20, which also shows that, when increasing a, the 
elastic piece rises faster than the inelastic one (compare to Fig. 18) and rapidly dominates 
over the latter at all dipole sizes, in agreement with our previous estimates. 

After replacing AP^m ~ APei in Eq. (4.47) and also using Eq. (4.48), it becomes clear 
that the diffractive probability coincides with the elastic probabihty, Eq. (4.46), which 
demonstrates Eq. (4.4) to the order of interest in ln(l//3). Although we have considered 
here only the simplest evolution of the projectile (namely, the formation of the qqg state), 
we expect the property (4.4) to hold also for smaller values of f3 (when higher Fock-space 
components come into play), for the reasons explained after Eq. (4.4). Assuming this to 
be the case, we conclude that, within the high-energy regime at ajo ^ 1 and Z <^ a^, the 




(4.64) 





(4.65) 



61 



10 F 



Tp- 1 i-ri— 1 r-q— 1 r-q— i r 



T 



1 r- 



Fluct. CT = 3 

Fluct. (7 = 5 

Fluct. (7 = 7- 



/' 



/ 




100 



Fig. 20. The contribution of the qqg component of the onium wavefunction to the elastic diffrac- 
tion probabihty as a function of r'^{Ql). The high-energy estimate (4.64) is shown for various 
values of a and compared to the mean-field prediction obtained with Eq. (4.8). 



diffractivc DIS cross-section is dominated by the elastic scattering of the onium, and thus 
is independent of (3 (for a given value of Y), within the limits of the leading-logarithmic 
approximation with respect to ln(l//3). To that accuracy, the corresponding cross-section 
has already been evaluated in Eq. (4.41). 

Let us conclude with a remark which applies to both the inclusive and the diffractive 
cross-sections and refers to their behaviour when increasing Y in the high-energy regime: 
At the end of Sect. 3, we have mentioned the 'rigidity' of the dipole amplitudes at high- 
energy, namely, their property to increase only slowly when increasing Y at fixed, but 
large Q^, in the weak scattering regime. Not surprisingly, the same 'rigidity' transmits to 
the DIS cross-section in the diffusive-scaling window defined by Eq. (4.17). For instance, 
the y-derivative of the diffractive cross-section in Eq. (4.41) is estimated as 



that is, it involves an additional suppression factor Z/cr^ -C 1 as compared to the usual 
behaviour expected in the context of the BFKL evolution. 




(4.66) 



which for Q'^ within the range indicated in Eq. (4.17) behaves like 




when CT < Z < (7^ (4.67) 
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5 Conclusions 



In this paper, we have formulated the inclusive and diffractive deep inelastic scattering 
of a virtual photon off a hadronic target in the high-energy limit and in the large-iV^ 
approximation. We have worked in a Lorentz frame in which the high-energy limit is 
achieved by boosting the hadronic target, which then evolves to very high gluon density, 
whereas the rapidity of the projectile (the virtual photon) remains finite and relatively 
small. Our main emphasis has been on a novel type of universal behaviour which emerges 
at sufficiently high energy as a consequence of saturation and gluon-number fluctuations 
in the target wavefunction. The essential feature of this new regime is the fact that the 
DIS cross-sections at fixed impact parameter are dominated by the physics of saturation 
up to very large values of Q^, well above the (average) saturation momentum (Qg) of the 
target (which itself increases exponentially with Y). And the distinguished signature of 
this behaviour is a new, dijfusive, scaling law, which should be obeyed by all the DIS 
cross-sections at sufficiently high energies, and which reflects the Gaussian nature of the 
probability distribution for the logarithm of the saturation momentum in the target. 

To establish this behaviour, we have relied on two main types of ingredients: (i) a set 
of factorization formulae which relate the DIS cross-sections to dipole-target scattering 
amplitudes, and (ii) the high-energy estimates for the latter, as determined by the 
Pomeron loop equations via the correspondence with statistical physics. 

Concerning point (i), let us recall here that we have restricted ourselves to diffractive 
processes in which the hadron undergoes elastic scattering, and to leading-logarithmic 
approximations for the high-energy evolutions of both the target and the projectile. Be- 
sides, we have assumed the projectile to remain dilute, and thus to obey BFKL evolution. 
Under these assumptions, we have shown that a physically transparent description of the 
diffractive processes, a la Good and Walker, emerges when viewing the collision in the 
frame in which the rapidity of the target is equal to the minimal rapidity gap. In this 
frame, the final state interactions automatically cancel out, and the projectile can be rep- 
resented as a distribution of dipoles which elastically scatter off the target. The effect of 
the collision is to destroy the coherence of the dipole superposition and thus transform 
the incoming photon into a collection of partons in the final state. 

Concerning point (ii), we should emphasize that the asymptotic functional form of the 
dipole amplitudes at high energy follows from very general and elementary considerations 
— the average amplitudes are dominated by saturated gluon configurations, for which 
the single-event amplitudes have reached their unitarity limit T = 1 — , and as such 
it should be very robust: It is likely that this dominant behaviour will not be modified 
by higher order perturbative corrections, nor by corrections in l/N^ (but this remains to 
be investigated). On the other hand, the energy dependencies of the average saturation 
momentum (Ql) and of the front dispersion cr^ are sensitive to the details of the evolution, 
and arc presently unknown for realistic values of as- This is why, in this paper, we have 
systematically measured in terms of {Q^, and Y in terms of cr^. 

Perhaps the most interesting, qualitative, conclusion of our analysis is the fact that, at 
sufficiently high energies and as a result of fiuctuations, the physics of gluon saturation 
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starts to be determinant at momenta Q"^ well above the saturation momentum (Q^), within 
a logarithmic distance ln(Q^/((5s)) ~ o"^ which increases linearly with Y . From previous 
studies of the mean field approximation, we already knew that saturation can manifest 
itself at momenta larger than Qs, via the phenomenon of geometric scaling. In the present 
analysis, we have shown that this phenomenon is preserved by fluctuations so long as the 
front dispersion is not too large (i.e., for o"^ ^ 1). What is however striking about the 
high-energy regime at o"^ > 1 is that the physics above (Q^) is not only influenced^ but 
actually dominated by rare configurations at saturation. This has tremendous physical 
consequences, for both the dipole amplitudes and the DIS cross-sections : 

• The strong correlation property {T'^{r)) ~ (T(r)), which is natural at saturation where 
(T) = 1, appears to be satisfied up to very large Q'^ = where (T) <C 1. 

• The amplitudes and the DIS cross-sections are very rigid, in the sense that they rise 
unusually slowly when increasing Y at fixed within the weak scattering regime. 

• At relatively large virtualities ^ (Qs)i cf. Eq. (4.17), the convolutions with the 
photon wavefunction yielding the DIS cross-sections are dominated by small sizes r ~ 
l/Q for the qq pair, as they select those gluon configurations in the target wavefunction 
which are at saturation on the resolution scale set by the virtual photon. 

• Similarly, the convolutions with the dipole probabilities yielding the average over the 
BFKL (onium) wavefunction of the projectile are dominated by dipole sizes of the order 
of their parent dipole — the original qq pair — and hence of order 1/Q. 

• In the high-energy limit, the DIS diffractive cross-section becomes purely elastic (at 
the level of the qq pair), the inelastic contribution being parametrically suppressed. 

• Within the high-energy regime defined by Eq. (3.11), all the amplitudes or cross- 
sections show diffusive scaling, that it, they depend upon Q"^ and Y only via the variable 
z/a with^ = ln(QV(g2)). 

• Within the diffusive scaling region alluded to above, the 'twist expansion' does not 
apply anymore: At large ^ (Qs)> cross-sections decay as a Gaussian of z/a, 
and not as a power of l/Q^. 

These properties are in sharp contrast with the corresponding predictions of the MFA (that 
we have also worked out in this paper, at least qualitatively, for the sake of comparison), 
and demonstrate the breakdown of the latter at sufficiently high energy. 

Note that the above conclusions depend upon the functional form of the dipole amplitudes, 
but not upon the precise values of coefficients A and Dfj.. But, of course, the latter are 
necessary in order to translate our results into physical units for Y and Q^. The small-x 
phenomenology at HERA suggests A ~ 0.2 0.3 [40,49], but we refrain ourselves from 
suggesting a similar, experimental, measure of the diffusion coefficient D^, as we do not 
beheve that the high-energy regime that we discuss here has been approached at HERA. 
(For instance, the diffractive data at HERA indicate the dominance of inelastic diffraction 
at small (3 [50], at variance with the high-energy behaviour that we predict here.) The 
conceptually proper way to compute A and Df^ for realistic values of (and within 
the present approximations: high energy, leading-order perturbative QCD, and large Nc) 
would be by solving the Pomeron loop equations of Ref. [14]. This would also allow one 
to explicitly test our present, asymptotic, predictions, and to study the preasymptotic 
behaviour, like the dependence upon the initial conditions or upon the impact parameter. 
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Such more detailed calculations should allow one to determine the physical region in 
which this high-energy behaviour should start to manifest itself and, especially, whether 
this region might be accessible at the LHC. 
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A The virtual photon wavefunction 

In the dipole frame we have adopted throughout this paper, the DIS cross-sections are 
factorized as the product of the photon wavefunction (7* — > qq) with the quark-antiquark 
interaction amplitude. In this appendix, we give the qq dissociation probabilities for the 
virtual photon. This well-known process can be computed in perturbative QED (to lowest 
order in a^^) and reads 

r; g) r = E ej { {v' + (1 - vr)QjKl{Qjr) + mjK^(g^r)} , (A.l) 

mv,r;Q)\'^^^e} {4gV(l - ^)^K^(g,r)} , (A.2) 

where g^ = ^(1 — t')g^ -|- rrij, ruf and e/ is the mass and the electric charge of the quark 
with flavor /, and Kq and Ki are the modifled Bessel functions. 

In order to get estimates of the DIS cross-sections, i. e. to estimate the convolution of this 
wavefunction with the qq interaction amplitude, we shall ignore the quark masses (these 
are important only for relatively small ~ mj) and use simple approximations for the 
Bessel functions, as quite common in the literature (see, e.g., [40] for similar manipula- 
tions). Specifically, since Kiy{x) decreases exponentially at large x, one can restrict the 
integrations in Eqs. (4.1)-(4.2) to values v and r such that Qfr < 1. Besides, for x ^ 1, 
Ki{x) ~ 1/x and Ko(x) ~ ln(l/a;), so the Bessel functions can be approximated as 

Ko(x) ~ 0(1 -x), Ki{x) r^-e{l-x), (A.3) 

where we ignore the overall normalization as well as the logarithmic singularity of Kq as 
X — > (the latter is innocuous within the relevant convolutions). Another simplification 
arises when performing the integral over v : the constraint v{l — v)Q'^ < 1 together with 
the fact that v{l — v) < l/A for 0<i'<l makes it natural to distinguish between 

(a) symmetric configurations, for which g^r^ < 4 ("small dipoles"), and 
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(b) a%ned je^ configurations, for which 4 < Q^r ("large dipoles"). 

In the first case, there is no restriction on the v integral, which is then dominated by 
symmetric values v ~ 1/2. In the second case, the dominant contributions correspond 
to the situation in which one of the two dipoles carries most of the total longitudinal 
momentum, that is, v is either close to zero, or close to one, which allows us to further 
simplify the respective integrand. As a result of such simplifications, the convolutions 
involving the photon wavefunctions can be finally estimated as shown in Eqs. (4.9)-(4.10). 



B The diffractive cross-section near j3 = 1 

In this Appendix we present the result for the differential diffractive cross-section per 
unit rapidity for /? ~ 1 as obtained when relaxing the leading-logarithmic approximation 
in ln(l//5). As compared to the formulae in the main text, this quantity represents the 
integrand which would give the elastic piece in the r.h.s. of Eqs. (2.1) if that quantity was 
computed beyond the leading-log approximation w.r.t. ln(l//?). Notice that, within the 
limits of that approximation, the elastic scattering implies (3=1, hence the corresponding 
differential distribution is just a 5-function ~ 5(/? — 1). But this distribution gets smeared 
after more accurately taking into account the kinematics for the diffractive production of 
the qq pair into which the virtual photon has fiuctuated. 

In the leading-logarithmic approximation with respect to ln(l/a;), one finds [36,37] 



^ J d^rd^r'dvv{l-v)Q{K}) e"^fir'-r) 



d%d\n{l/(3) At: (5 ^ 

E $^(^,r,Omr,6))^,,,(r(r',6))y,,, (B.l) 

a=L,T 

where lij = v{l — v)M'j^ — mj. The photon wavefunctions and are given by 

4(^, r, r') = ^ej((^^ + (1 - ^ K,(Q,r) Ki(Q/r') 

+ mjKo(g/r)Ko(g/r')) (B.2) 
(we use the same notations as in Appendix A) and respectively 

^Uv, r, r') = ^ e] AQ\\l - vf MQfv) Ko(Q/r') (B.3) 
with r = |r| and r' = \r'\. Note that these functions are such that 

Y^^i{v,r,r) = \ra{v,r-Q)\^ (B.4) 
/ 

is the photon wavefunction squared discussed in Appendix A. Note furthermore that the 
r.h.s. of (B.l) depends upon ^, via h.). (Recall that (3 = Q'^/{Q'^ + Af|).) 
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The cross-section (da^fg/d^fe) is obtained by integrating (B.l) from ln(l//5) = up to 
Y — Fg™™. If this upper hmit is relatively large, one recovers the qq part of our expression 
in the main text, cf. Eq. (4.6) (with Y^^^ renoted as Vgap, for simplicity) : 



'iY,Y,^„Q') = fd'rfdv Y: mv,r;Q)mT{r,b)ry^^^ . (B.5) 

a=L,T 

This can be seen by changing the integration variable to K.'j : the typical values of r 
and r' are set by the photon wavef unctions, and thus are of order l/Qf. On the other 
hand, the integral over fixes the difference r — r' to be of order with kJ^"^ 

determined by the maximal diffractive mass Mx (the one corresponding to the minimal 
rapidity gap l^gap)- When k,J^^ ^ Qf (in practice, this condition amounts to Af| ^ Q^), 
the difference r — r' is small as compared to both r and r', so one can approximate the 
integral / d^Kj e*'*/'^'' "'^^ by S^'^\r' — r), which then yields (B.5). 

Note finally that the differential cross-section is proportional to the diffractive structure 
function F^'^{x,Q^, P), which is measured at HERA: 
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